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SUMMARY 


Structural  response  predictions  have  been  made  for  two  Inportant  areas  of 
a short-takeoff-and-landing  (STOL)  aircraft.  The  method  was  developed  to 
significantly  Improve  environmental  predictions  compared  to  those  used  In 
the  past.  A mathematically  rigorous  spectral  analysis  approach  was 
developed  that  simulated  the  structure  with  a finite-element  model  which 
used  correlated  and  calculated  acoustic  Input  data  for  the  forcing  function. 

The  structural  vibration  analyses  were  successful  1n  predicting  operating 
levels  and  describing  the  frequency  content  of  responses  at  selected 
locations  on  the  structure.  Comparisons  of  predicted  and  measured  data 
show  that  the  method  developed  and  described  here  may  be  used  to  predict 
complex  structural  response  to  jet  engine  excitation. 

A method  was  also  developed  for  the  prediction  of  the  external  acoustic 
environment  of  STOL  aircraft  with  upper-surface-blown  (USB)  flaps.  The 
method  Is  described  In  detail,  and  comparisons  are  given  between  predicted 
and  actual  measurements.  The  method  gives  good  results  and  represents  a 
significant  Improvement  over  previous  acoustic  prediction  methods. 


SECTION  1 
INTRODUCTION 


1.1  Program  Objectives 

The  first  objective  of  the  program  was  to  use  available  data  of  the  vibration  and  acoustic 
characteristics  on  Short  Takeoff  and  Landing  (STOL)  aircraft  to  predict  the  resulting 
aircraft  structural  vibration  levels. 

A second  objective  was  to  develop  a method  of  predicting  the  external  acoustic  levels 
for  a STOL  aircraft  and  to  use  these  predictions  as  Inputs  to  the  structural  vibration 
analysis  program. 

The  two  areas  chosen  for  detailed  study  of  both  objectives  were  the  wing/ flap  structure 
and  the  fuselage  section  adjacent  to  the  wing  root,  upper  surface. 

Use  of  the  methods  developed  will  provide  environmental  vibration  predictions  in  all 
areas  of  STOL  type  aircraft. 

Additionally,  parametric  studies  were  made  of  STOL  aircraft  from  50,000  lb  to  1,000,000 
lbs  for  structural  vibration  prediction  levels.  These  values  were  compared  to  the 
vibration  criteria  of  MIL-STD  810  C, 

IJZ  Program  Definition 

Phase  1 of  the  STOL  Program  was  divided  into  10  tasks.  These  includei 


Task  1 
Task  2 
Task  3 
Task<f 
Task  5 
Task  6 
Task  7 
Task  8 


Program  Definition 

Method  Development  - Flap  Structures 

Apply  YC-14  USB  Flap  Data 

Comparison  to  Flap  Data  Tests 

Discrepancies  and  Refinement  of  Flap  Prediction 

Fuselage  Structure  Method  Development 

Apply  YC-14  Fuselage  Data 

Comparison  to  Fuselage  Data  Tests 


Task  9 
Task  10 


Discrepancies  and  Relinement  of  Fuselage  Predictions 
Acoustic  Field  Prediction  Development 


Phase  II  was  divided  into  four  tasks.  These  include: 

Task  11  Parametric  Studies  of  Vibration  Response 

Task  12  Noise  Field  Parametric  Predictions 

Task  13  Compare  Predictions  to  Estimated  Test  Specs 

Tsisk  14  Report  Preparation 
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PHASE  I - DEVELOP  PREDICTION  METHOD 


SECTION  II 

TECHNICAL  ANALYSIS  METHOD  DEVELOPMENT 

The  technical  analysis  for  Phase  I consisted  of  three  parts:  1)  eigenvalue  analysis  of  a 
finite  element  structural  model,  2)  definition  of  the  acoustic  environment,  and  3) 
random  harmonic  analysis.  The  data  flow  for  the  analysis  appears  in  Figure  1, 

For  the  finite  element  models  used  in  the  analysis,  two  USB  flap  models  and  three 
fuselage  models  were  developed.  These  are  described  in  greater  detail  in  Sections  III  and 
IV.  Inputs  to  the  finite  element  analysis  included  model  geometry  and  degrees  of 
freedom  of  an  array  of  nodes;  definition  of  a system  of  structural  elements,  fixity  of 
each  element,  section  properties,  material  properties  and  assumed  structural  damping 
values.  The  Structural  Analysis  Program,  SAP  IV,  (Reference  1),  written  lor  the  CDC 
6600  computer,  was  used  for  the  finite  element  model  eigenvalue  analysis,  Outputs 
included  modal  frequencies  and  modal  displacements  at  each  node  location.  The  acoustic 
environment  used  for  the  analysis  in  Phase  I was  based  on  available  test  data  which 
consisted  of  power  spectral  density  plots  of  sound  pressure  vs  frequency  at  several 
microphone  locations.  A scheme  was  devised  for  interpolating  and  extrapolating  the 
data  to  each  panel.  This  interpolation  scheme  is  discussed  in  Section  III.  The 
interpolation  was  accomplished  through  a simple  computer  program  which  generated 
output  data  in  the  complex  matrix  format  required  for  subsequent  analyses.  Real  and 
imaginary  parts  of  each  matrix  element  were  generated,  with  one  matrix  being 
generated  for  each  input  frequency.  The  matrix  size  is  equal  to  the  number  of  panels. 
The  diagonal  elements  (real)  are  values  of  power  spectral  density  at  each  panel.  The  off- 
diagonal  terms,  representing  cross-spectral  density  terms,  were  included  in  the  first  set 
of  calculations  for  the  USB  flap  response  and  compared  to  the  results  from  calculations 
where  the  off-diagonal  terms  were  set  equal  to  zero.  The  results  differed  by  only  10%, 
so  subsequent  computations  were  made  using  only  the  diagonal  elements. 

For  the  structural  response  analysis,  it  was  necessary  to  know  the  modal  displacements 
at  the  centroids  of  each  panel.  A simple  computer  program  was  written  to  interpolate 
the  mode  shapes  from  the  finite  element  analysis.  This  interpolation  consisted  of  simple 
arithmetic  averaging  of  the  modal  displacement  at  each  of  the  four  corners  of  a panel. 


The  modal  data  and  the  acouitic  envlronmant  data,  were  Input  to  the  Random  Harmonic 
Analyais  Program  (TEV  156,  Reference  2).  Also  Input  were  values  of  generalized 
stlffnessi  structural  damping,  modal  displacements  at  output  stations  (corresponding  to 
accelerometer  locations  for  which  test  data  are  available),  and  a list  of  frequendas  at 
which  output  data  w«e  desired.  Details  of  the  solution  technique  employed  appear  in 
Reference  2.  Briefly,  the  program  solves  Equation  1 of  Reference  2 which,  when 
modeled  for  this  application  la* 


where 


CmJ 

CmJ 

[C,] 

Cf] 


Generalized  stiffness  matrix  (size  m x ' 
Generalized  inertia  matrix  (size  m x m. 
[♦]^[a]  ■Fordng  function  matrix  (size  m x n) 
Matrix  of  modal  dlspaicemcnta  (size  n xm) 
Matrix  of  panel  areas  (size  n x n) 

Matrix  of  generalized  coordinates  (size  m x 1) 
Matrix  of  generalized  accelerations  (size  m x 1) 
Structural  damping  coeffident 

matrix  of  panel  pressures  (size  n x i) 

Number  of  modes 
Number  of  panels 


The  load  equations  follow  the  same  format  as  the  equations  of  motion. 

where  | | ■ Matrix  of  accelerations  at  output  stations  (size  k x 1) 

k ■ Number  of  output  stations 


C « (l/g)  ■ Matrix  of  coeffldents  (size  K x m) 

I > Gravitational  acceleration 

[^a]  ” Matrix  of  modal  displacements  at  output  sta 


Matrix  of  modal  displacements  at  output  stations  (size  k x m) 
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The  program  performs  Laplace  transformations  on  Equations  1 and  2,  resulting  in, 


[ [“l]  20*S^[“3]]{5}=[c3]f  (p)} 

{uz)}  = S^[!r3]{^} 


where:  s 

{5t 

|UP)1 


Laplace  operator 

Laplace  transformed  coordinate  matrix 
Laplace  transform  of  panel  pressures 
Laplace  transform  of  output  accelerations 


(3) 

(4) 


Generalized  coordinate  and  load  frequency  response  functions  are  obtained  by  solving 
Equations  3 and  4. 


The  program  then  employs  the  technique  of  generalized  harmonic  analysis,  viz., 


♦o  “ [tJ  [cpsd]  {t} 


(5) 


where: 

LtJ 

{t} 

[cpsd] 


Output  power  spectrum  at  a specific  station  at  frequency 

Row  matrix  of  output  frequency  response  to  a sinusoidal  force  of 

f K 

frequency  0)  acting  at  the  i excitation  point  (size  1 x n) 
Transpose  of  T 

Matrix  of  cross-power  spectral  densities  at  frequency  U)  (size  n x 
n) 


Equation  (5)  is  solved  for  each  frequency  of  interest  and  for  each  output  station.  The 
root-mean-square  response  is  obtained  by  the  expression: 


00 
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SECTION  m 


FLAP  STRUCTURE  VIBRATION  PREDICTION 

This  section  of  the  study  is  limited  to  the  YC-14  aircraft  USB  flap  structure,  since 
significant  amounts  of  ground  and  flight  test  data  were  available  for  comparisons. 

Both  vibration  and  acoustic  data  were  obtained  from  ground  and  flight  tests.  These  data 
were  recorded  simultaneously  for  several  conditions  to  enable  correlation  studies  to  be 
made  of  the  acoustic  input  and  the  structural  response.  The  flap  area  received  the 
highest  input  energy  on  the  aircraft  and  thus  determination  of  the  structural  response  in 
this  area  was  of  paramount  interest. 

3.1  USB  Flap  Model  Development 

The  first  Upper  Surface  Blowing  (USB)  flap  model  was  considered  as  two  separate  panels 
(plates);  the  main  USB  flap  and  the  aft  USB  flap. 

To  solve  the  frequency  determination  of  plates,  an  approximate  solution,  using  the 
Rayleigh  principle,  was  used  with  the  Warburton  Method  (Reference  3)  where  the 
coefficients  In  the  frequency  equation  were  given  for  several  different  boundary 
conditions. 

Calculations  were  made  for  a plate  simulating  the  YC-l<i  main  USB  flap.  The  initial 
boundary  condition  used  was  condition  15  on  page  375  of  Reference  3 where  the  forward 
edge  was  assumed  to  be  fixed. 

The  comparison  of  calculated  and  experimental  frequencies  indicated  a large  discrepancy 
between  use  of  Warburton's  prediction  for  a flap  aspect  ratio  of  1.716  and  1.225,  i.e., 
ratios  of  span  to  chord  dimertsion  of  main  and  aft  USB  flaps. 

In  reviewing  the  frequency  data  of  both  the  YC-14  airplane  main  and  aft  flaps,  it  was 
noted  that  the  frequency  of  both  were  nearly  identical.  This  fact  clearly  leads  us  to 
believe  the  two  flaps  act  very  nearly  as  a single  unit.  With  this  in  mind,  the  decision  was 
made  to  calculate  Warburton's  frequency  predictions  using  the  combined  flap  assembly  as 
a single  unit,  where  the  aspect  ratio  was  .2941. 
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The  frequency  comparisons  indicated  a good  first  approximation  but  not  sufficiently 
accurate  to  warrant  use  of  this  model  for  response  calculations  using  acoustic  excitation 
of  the  flap.  In  the  Interest  of  a more  rigorous  math  model  and  an  overall  calculation 
method,  the  decision  was  made  to  build  a finite  element  model.  With  such  a model,  the 
details  of  the  structure  could  be  more  accurately  described  and  response  of  given 
locations  compared  to  the  flight  test  data 

Two  finite  element  USB  flap  models  were  then  used  In  the  development.  Model  I was  a 
simple  plate  finite  element  model  and  was  found  to  lack  some  definition.  As  a result,  a 
more  detailed  model  was  formulated  in  Model  II  and  provided  more  detailed  results. 


3.2  Finite  Element  USB  Flap  Model  I-Definitlan 

Model  I was  fashioned  after  the  plate  model  that  was  studied  earlier  with  the  Warburton 
calculations.  This  model  was  formulated  as  drawn  in  Figure  2.  The  dimensions  for  the 
model  nodal  points  were  selected  as  shown  in  Figure  3. 

Model  I was  then  input  to  the  SAP  IV  program  with  the  input  format  as  given  in  the 
print-out  of  the  data  card  image  of  Figure  4.  The  print-out  of  the  frequencies  for  the 
first  10  modes  was  listed  in  Figure  5.  This  data  was  plotted  for  comparison  to  test  data 
in  Figure  6. 

Model  I was  formulated  after  the  mass  was  matched  to  the  actual  flap  weight  of  837 
pounds.  The  Warburton  calculations  had  the  dimensions. 


L > 

204  in. 

(length) 

W s 

60  in. 

(width) 

t.  = 

.94  in. 

(thickness) 

11 

Oa  Ib/in^ 

(Mass  density) 

386 

Thus,  W^ 

* Lw  t„ 

^0  s 1130.36  lbs. 

or, 

837 

m .7273,  was  the  correction  factor  for  mass. 

TTiro 

The  SAP  IV  run,  using  Warburtons  thickness  as  a first  Input,  gave  a frequency  ratio  from 

experiment  of  7.64  as  seen  in  the  experimental  data,  i.e., 
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These  values  of  thickness  and  mass  density  were  then  used  as  input  to  the  finite  element 
program,  SAP  IV.  The  resulting  model  frequencies  are  shown  In  Figure  3 for  the  first  ten 
modes.  Comparison  with  measured  values  is  shown  in  Figure  6. 


Modal  plots  have  been  generated  from  these  SAP  IV  runs,  using  the  tabulated  data  as 
shown  in  Figure  7.  The  actual  mode  shape  plots  were  obtained  from  this  tabulation  and 
are  shown  in  Figure  8 thru  13.  It  can  be  seen  that  the  modal  definitions  are  reasonable  fcr 
the  density  of  the  data  points  taken.  This  model  was  thus  felt  to  represent,  in  a realistic 
manner,  the  USB  flap  and  warranted  the  application  of  an  Input  excitation  for  harmonic 
analysis. 

3.3  Harmonic  Analysis  of  Model  1 

The  USB  flap  response  power  spectrum  for  surface  element  g due  to  the  pressure 
excitation  forces  on  ail  surface  elements  is  of  the  form: 
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«,(u) 

“ E E *!,<“>  Tj,  («•»  - 1.T.J  [cWd]  LtJ  (« 

where 

♦q(u) 

> Response  power  spectrum  for  element  q 

8 Cross-power  spectra  (CPSD)  of  the  excitation  forces  on 
elements  1 and ) 

V“’ 

8 Output  frequency  response  at  element  q to  a unit  lurce 

AW 

acting  at  the  1 element. 

The  excitation  points  were  determined  by  dividing  the  USB  fiap  upper  surface  into  i6 
paneis  of  equai  area  as  was  shown  in  Figure  with  a given  pressure  acting  over  each 
panei.  The  power  spectra  and  cross>power  spectra  for  the  pressures  acting  on  the  paneis 
were  extrapoiated  from  the  data  that  was  obtained  from  the  four  acoustic  sensors 
located  as  shown  in  Figure  14.  The  extrapolation  of  the  data  for  each  panel  for  which 
there  were  no  measured  data  was  accomplished  using  the  following  formuiast 

Measured  Data 

*10,  10}  *11,  llj  *13,  13j  *14, 14|  *10,  14j  *11,  14j  *13,  14 

Extrapolation  Formula 

*9,  9 * *10,  10 
12,  12  13,  13 

*13,  13  • *16,  16  “*14,  14 

*i,l  “ *(i+8Xi+8)  ' • • • * 

*l,(i+l)  * *13,14  , where  Ul, ...  7 and  1=9, . . .13 

*i  (i+2)  " *11  14  ’ 6 1=9, . . .14 


I 
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^l,(l+3)  “ ^11,14  * • • • 3 «nd  i«9, ...  13 

* ^10  1^  ^ wh#r®  IbXi  m m m ^ 40)ci  ...  12 

^i,a+7)  * *i3,i» » *■**  • • • * 

^l,(l-4)  ■ *13,14  » * • • * 

*l,(l+9)  • *13,14  » • • * ^ 

A map  of  the  uppw  triangl.  of  th.  CPSD  matrix  Is  shown  In  Figure  13.  The  lower  part  of 
the  triangle  Is  the  complex  conjugate  of  the  upper  triangle. 

The  ec^ations  of  motion  and  load  equations  were  based  on  modes  calculated  in  SAP  IV. 
Equations  of  Motion 

*WH*«c[m.](,)  .[c,][p) 

Generalized  stiffness  matrix  (size  m x m) 

Generalized  Inertia  matrix  (size  m x m) 

[^]  [ A j ■ Forcing  function  matrix  (size  m x n) 

Matrix  of  modal  dlspalcements  (size  n xm) 

Matrix  of  panel  ar  eas  (size  n x n) 

Matrix  of  generalized  coordinates  (size  m x 1) 

Matrix  of  generalized  accelerations  (size  m x 1) 

Structural  damping  coefficient 


[mJ, 

where 
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matrix  of  panel  preuures  (size  n x 1) 
Number  of  modes 
Number  of  panels 
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An  example  of  a portion  of  the  output  spectrum  values  for  the  three  accelerometer 

positions  are  shown  on  the  computer  printout  sheet  in  Figure  16.  These  results  have  been 

2 2 

tabulated  and  converted  from  G /RAD/sec  to  G /Hz  as  shown  in  Figure  17. 

The  solutions  were  for  YC-14  flight  condition  where  the  test  conditions  were;  altitude 
7,620  feet,  speed  216  ft/sec,  Nj  of  3,09S  RPM,  and  USB  flap  setting  of  40  degrees. 


The  resulting  RMS  accelerations  for  the  three  accelerometer  locations  were  as  follows, 
tabulated  below: 


ACCELEROMETER  NO. 

DENSE  CPSD 
MATRIX 
g = .03 

^^RMS^ 

DENSE  CPSD 
MATRIX 

8 = .01 

^^RMS^ 

DIAGONAL  CPSD 
MATRIX 

^^RMS> 

1421 

1.82 

3.02 

1.67 

1417 

2.93 

4.82 

2.53 

1428 

1.38 

2.38 

1.18 

LOAD  EQUATIONS 


s frequency  responses  for  acceleration 
a matrix  of  modal  values  at  accelerometer  locations 
= The  matrix  of  generalized  coordinates,  each  column  representing 
the  response  to  a unit  pressure  acting  on  one  of  the  excitation  panels. 


[^] 

Ct] 

[♦o] 

c^] 

Procedure 


At  each  frequency,  the  dynamic  analysis  computer  program  solves  for  the  response, 
^tJ  , and  performs  the  calculations  as  shown  in  Equation  (8).  Each  output  spectrum  is 
integrated  over  the  specified  frequency  range  (26  to  193  Hz)  to  obtain  the  RMS  value  of 
the  accelerometer  response.  The  cross-power  spectra  were  enriched  by  extrapolating 
the  measured  data  as  described  in  the  list  of  extrapolation  formulas  given  previously. 
The  basic  equations  of  motion  using  the  generated  modal  values  from  the  SAP  IV  finite 
element  program  yielded  the  generalized  coordinates  of  the  response  to  a unit  pressure 


II 


iH^ 


) 

i 
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acting  on  the  excitation  panels.  Once  the  response  coordinates  were  obtained  the 
frequency  responses  could  be  determined  as  indicated  in  Equation  (8). 

3.4  Comparison  of  Model  I and  Test  Data 

The  following  experimental  values  were  obtained  for  the  frequency  range  used  in  the 
calculations  (26-193  Hz): 


Acceleration  No. 

‘^RMS 

1421 

1.39 

1417 

1.20 

1428 

0.40 

Several  points  must  be  discussed  before  any  general  conclusions  can  be  drawn.  Model  1 
does  not  contain  the  hard  points  of  flap  actuation  attachment,  where  the  accelerometers 
were  located.  Neither  does  this  math  model  simulate  the  heavy  spars  In  the  flaps  where 
those  attachment  points  are  located.  From  this  lack  of  simulated  tie-down  or  hard 
points,  we  would  expect  the  calculated  results  to  lack  definition. 

The  initial  comparisons  for  the  USB  flap  setting  of  40°  were  for  three  accelerometer 
locations  where  calculations  for  g = .03  were  used  with  the  "dense"  cross-power  spectral 
density  (CPSD)  matrix.  Calculations  for  these  accelerometers  using  the  diagonal  CPSD 
matrix  gave  results  that  were  within  10%  of  the  results  from  the  "dense"  CPSD  matrix. 
Thus,  only  the  diagonal  CPSD  matrix  was  used  for  subsequent  comparisons. 

The  response  values  for  the  YC-14  USB  flap  were  calculated  for  additional  damping 
values  of  .06,  .09,  .12  and  .13|  see  Figure  18.  The  total  damping  value  for  the  USB  flap 
of  0.13  represents  the  best  match.  This  value  would  Include  structural  damping  as  well 
as  aerodynamic  damping. 

The  response  values  for  the  three  USB  flap  accelerometers  for  an  assumed  structural 
damping  value  of  g = .13  are  shown  In  Figures  19,  20,  and  21. 
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3.5  Finite  Element  USB  Flap  Model  n-Definition 

The  results  from  the  Model  I were  encouraging  but  did  Indicate  a better  model  was 
• needed  for  more  detailed  response  predictions.  Thus,  Model  II  was  formulated  and  much 
of  the  actual  USB  flap  structure  was  simulated  including  the  attachment  points  of  the 
two  hydraulic  actuators  at  node  points  40-73  and  45-74  of  Figure  22.  Details  of  this 
model  appear  in  Figures  23  through  25.  The  refinements  Incorporated  into  Model  II 
which  did  not  exist  in  Model  I,  include  the  followings  (I)  The  number  of  nodes  has  been 
increased  by  a factor  of  nearly  3f  (2)  Actual  geometry  is  better  represented;  (3)  .Spars, 
- ribs,  leading  edges,  and  trailing  edges  are  represented  as  beam  elements;  (4)  Hinges  are 
free  to  rotate;  (5)  Actuators  are  simulated  by  truss  elements  (6)  Actual  material 
properties  are  used  In  the  simulation  of  the  various  structural  components. 

The  finite  element  model  Includes  78  nodes,  74  of  which  are  on  the  flap  and  4 of  which 
represent  actuator  attachment  points.  Structural  elements  include  64  beams,  59  plates, 
and  two  truss  elements.  Although  the  skin  thickness  and  material  vary  over  the  flap,  a 
constant  equivalent  aluminum  plate  thickness  was  used  In  the  model.  Since  the  flap  Is 
actually  a three-dimensional  structure,  it  was  necessary  to  use  equivalent  plate  thickness 
as  described  previously  for  Model  I.  Using  the  SAP  IV  program,  following  an  initial 
eigenvalue  analysis  and  mode  shape  inspection,  a procedure  similar  to  that  described 
previously  was  employed  to  obtain  a frequency  match. 

The  list  of  the  first  13  natural  modes  of  the  USB  flap  Is  given  in  Figure  26.  These  mode 
shapes  resulting  from  the  finite  element  analysis  have  been  plotted  in  Figures  27  thru  39 
and  illustrate  the  complexity  of  the  USB  flap  vibrational  response  thru  a frequency  of 
approximately  300  Hx. 

3X  Harmonic  Analysis  of  Model  n (78  Node  ModeO 

The  harmonic  analysis  of  the  USB  flap.  Model  II  was  made  using  the  acoustic  excitation 
measured  by  the  four  acoustic  microphones,  M35,  M37,  M40  and  M41,  as  shown  in  Figure 
40  as  the  input  forcing  function.  This  acoustic  data  for  the  four  microphones  is  given  in 
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Figures  through  44.  These  data  were  obtained  from  the  YC-14  in  a STOL  condition, 
altitude  7620  feet,  speed  216  ft/sec  and  a USB  flap  angle  of  40°.  The  diagonal  CPSD  was 
used  for  the  Input  excitation  and  the  data  interpolated  to  cover  the  entire  model.  The 
Random  Harmonic  Analysis  was  then  used  with  thr^  USB  flap  finite  element  Model  II 
simulation  to  determine  the  response  values  for  the  three  locations  corresponding  to 
accelerometers  No.  1417,  1421  and  142S  shown  with  asterisks  in  Figure  43.  The  resultant 
response  predictions  have  been  tabulated  in  Figure  46  for  three  different  assumed 
damping  values.  A sample  page  from  the  computer  print-out  for  the  accelerometer 
response  prediction  with  g s .06  Is  given  in  Figure  47. 

3.7  Comparison  of  Predictions  for  Model  11  With  Test  Data 

The  predictions  for  G = .13  for  all  three  accelerometer  locations  have  been  plotted  In 
Figures  48,  49  and  30.  The  actual  test  data  from  three  accelerometers  that  were  located 
on  hard  structure  where  actuators  were  attached  are  also  shown.  The  detailed  structure 
at  these  locations  was  not  completely  simulated  but  the  results  show  levels  that  were 
very  representative  of  the  high  environment  associated  with  this  area.  The  frequency 
content  also  Is  noted  to  be  indicative  of  the  frequency  range  in  the  higher  environment. 
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SECTION  IV 


FUSELAGE  STRUCTURE  VIBRATION  PREDICITON 
4.1  Fuselage  Finite  Eiement  Model  Development 

The  fuselage  area  of  interest  Is  shown  in  Figure  51.  The  primary  considerations  in  the 
finite  element  structural  modeling  werej  (1)  The  need  to  cover  a broad  frequency  range 
(25  - 1000  cps);  (2)  Computer  resource  limitations;  (3)  Computation  costs.  It  was 
determined  that  a single  finite  element  model  would  not  be  adequate  for  the  entire 
frequency  range.  Three  models  were  developed,  one  for  the  low  frequency  range  (25  to 
100  cps),  one  for  the  intermediate  range  (100  to  200  cps),  and  one  for  the  high  frequency 
range  (above  200  cps).  Figures  52  through  54  show  the  nodal  grids  of  each  model  relative 
to  the  actual  fuselage  structure. 

Details  of  each  model  appear  in  Figures  55  through  67.  For  the  low-frequency-range 
model,  nodes  were  located  at  the  Intersections  of  every  third  frame  and  every  fourth 
stringer  as  shown  In  Figure  52. 

A finer  nodal  grid  was  selected  for  the  Intermediate-frequency-range  model  as  shown  in 
Figure  53.  Nodes  were  located  at  the  intersections  of  each  frame  and  stringer. 
Additional  nodes  were  located  on  each  stringer  at  points  midway  between  frames.  The 
nodal  density  was  24  times  that  of  the  low-frequency-range  model. 

The  high-frequency-range  model  was  represented  by  a nodal  grid  as  shown  in  Figure  54. 
The  central  portion  of  the  model  employed  a finer  grid  than  the  outer  portion  to  give 
better  definition  in  area  of  measurements.  Nodes  were  placed  at  each  frame/strlnger 
Intersection.  Additionally,  in  the  streamwise  direction,  seven  rows  of  nodes  were  placed 
equally  spaced  between  successive  frames.  Nodes  were  placed  midway  between  stringers 
in  the  outer  portions  of  the  model  and  three  equally  spaced  rows  of  nodes  between 
stringers  in  the  central  portion.  In  the  outer  portions,  the  grid  density  was  S times  that 
of  the  Intermediate-frequency-range  model  and  192  times  that  of  the  low-frequency- 
range  model.  In  the  central  portion,  the  grid  density  was  16  times  that  of  the 
intermediate  model  and  384  times  that  of  the  low-frequency-range  model,  giving  much 
increased  definition  for  determination  of  the  higher  frequency  modes. 
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Nodes,  coordinates,  and  structural  elements  for  the  low  frequency  range  Model  I are 
shown  In  Figures  55  through  58.  The  complete  model  simulates  a half-cylinder  section  of 
the  fuselage  spanning  16  frames.  Employing  the  grid  described  previously,  the  model 
consists  of  66  nodes.  Between  each  successive  set  of  stringer  nodes  is  a beam  element 
with  four  times  the  cross-sectional  area  and  four  times  the  moment  of  inertia  of  a 
stringer.  Between  each  successive  set  of  nodes  in  a tangential  direction  is  a beam 
element  with  three  times  the  section  properties  of  each  frame.  The  total  number  of 
beam  elements  is  116.  Plate  elements,  55  in  number,  are  located  between  each  set  of  4 
adjacent  nodes.  Element  material  and  section  properties  appear  in  Figure  67>  With  the 
exception  of  nodes  1-6  and  nodes  61-66,  each  node  was  given  two  translational  degrees  of 
freedom,  Y and  Z,  and  three  rotational  degrees  of  freedom.  Symmetrical  boundary 
conditions  were  imposed  on  nodes  1-6  and  nodes  61-66,  i.e.,  these  nodes  were  constrained 
from  displacement  in  the  Z direction  and  from  rotation  about  the  X and  Y axes.  The 
corner  nodes,  1,  6,  61,  and  66  were  constrained  from  any  translational  or  rotational 
motion. 

Nodes,  coordinates,  and  structural  elements  for  the  intermediate  frequency  range  Model 
11  are  shown  in  Figures  59  through  62.  The  model  simulates  a section  of  fuselage  defined 
by  four  frames  and  five  stringers.  The  model  Includes  35  nodes.  Beam  elements  with 
section  properties  equal  to  those  of  the  actual  structure  are  located  at  the  stringer  and 
frame  locations.  Plate  elements  are  also  included.  Material  and  section  properties 
appear  in  Figure  67.  Each  node  was  given  two  rotational  degrees  of  freedom,  about  the 
X and  Y axes.  All  but  the  corner  nodes,  1,  7,  29,  and  35,  were  given  a Z translational 
degree  of  freedom. 

Details  of  the  high  frequency  range  Model  ill  are  given  in  Figures  63  through  66.  The 
model  Includes  81  nodes,  48  beams,  and  64  plates.  The  degrees  of  freedom  were  the 
same  as  for  the  intermediate  frequency  range  model,  i.e.,  Z translation  and  X and  Y 
rotation  for  all  but  the  corner  nodes  (1,  9,  73,  and  81)  which  were  constrained  from  Z 
translation. 

The  above  descriptions  have  been  given  to  Indicate  the  degree  of  detail  that  would  be 
used  to  cover  given  frequency  ranges. 


4.2  Low  Frequency  Fuselage  Model  1 


The  SAP  IV  finite  element  program  was  used  to  predict  the  lower  frequencies  of  Interest 
of  the  Fuselage  Model  I.  A tabulation  of  the  first  20  frequencies  has  been  listed  In 
Figure  68.  The  mode  shape  of  the  center  portion  of  the  model  have  been  plotted  and 
serve  to  Indicate  the  response  of  such  a structure  nearly  free  of  the  end  constraints.  In 
this  case,  the  model  was  supported  with  the  four  extreme  corners  clamped.  These  node 
points  were  1,  6,  61  and  66  as  seen  In  Figure  35.  The  mode  shapes  of  the  center  portion, 
defined  by  node  points  19,  24,43,  and  48  have  been  plotted  in  Figures  69,  70,  7 1 and  72. 

The  acoustic  excitation  for  this  portion  of  the  fuselage  was  available  only  at  a limited 
number  of  transducer  locations.  Microphones  M6,  M13,  M16,  M18  and  M20  were  used  and 
Interpolation  and  extrapolation  were  necessary  so  as  to  cover  the  entire  model.  The 
locations  of  the  microphones  have  been  shown  In  Figure  73.  The  condition  chosen  for  the 
response  study  was  for  maximum  engine  thrust  during  ground  run-up.  The  acoustic  data 
was  plotted  in  Figures  74,  75,  and  76.  The  summary  of  the  acoustic  data  was  listed  In 
Figure  77  with  extrapolated  data  for  two  microphones  M18  and  M20. 

The  acoustic  PSD  data  was  then  used  to  excite  the  low  frequency  fuselage  model  of 
Figure  52  and  the  response  from  the  Harmonic  Dynamic  Analysis  Program  was  obtained 
at  locations  on  the  finite  element  model  corresponding  to  locations  where  actual 
accelerometers  were  located  in  the  YC-14  airplane.  The  accelerometers  were  located  as 
shown  In  Figure  73. 

The  calculated  responses  of  the  low-frequency  fuselage  model  have  been  plotted  at 
corresponding  locations  on  the  model  for  the  accelerometer  locations  of  A58  (stringer), 
A59  (body  frame)  and  A61  (skin,  center  of  panel)  and  compared  to  the  PSD  data  from  the 
actual  accelerometer  responses.  These  comparisons  have  been  plotted  In  Figures  78,  79, 
and  80. 


4 J Mid  - Frequency  Fuselage  Model  n 


The  mid-frequency  fuselage  Model  11  was  shown  In  Figure  53.  The  eigenvalues  for  the 
first  8 modes  are  shown  in  Figure  81  and  cover  a frequency  range  from  60  to  350  hz.  The 
model  was  pinned  at  the  four  corners  at  node  points  1,  7,  29  and  35.  The  mode  shapes 
have  been  plotted  in  Figures  82  thru  89.  The  acoustic  excitation  of  the  model  was  again 
obtained  from  extrapolated  data  from  microphones  M6,  M13,  M16,  and  M20.  The  same 
Flight  Test  Condition  was  used  to  calculate  model  response.  The  orientation  of  the  Mid- 
Frequency  Fuselage  Model  II  in  relation  to  the  low-frequency  fuselage  Model  I is  shown 
in  Figure  90.  In  addition,  the  location  of  the  microphones  and  response  measuring 
accelerometers  are  shown  in  this  same  Figure  90. 

The  response  of  the  fuselage  local  section  Model  II  has  been  plotted  also  in  Figures  78, 
79,  and  80  lor  the  three  accelerometer  locations. 

4.4  High  Frequency  Fuselage  Model  III 

The  High-Frequency  Fuselage  Model  III  was  shown  in  Figure  54.  The  detail  node 
positions  have  been  defined  in  Figure  63,  the  coordinates  ip  Figure  64  , the  beam 
elements  in  Figure  65  and  the  plate  elements  in  Figure  66.  The  SAP  IV  program  gave  the 
first  20  frequencies  from  269  to  987  Hz  as  listed  in  Figure  91.  The  mode  shapjss  have 
been  plotted  in  Figures  92  thru  95. 

The  acoustic  excitation  was  the  same  as  for  the  previous  models  with  the  Harmonic 
Analysis  results  of  the  response  predictions  for  the  stringer  (A58)  and  the  skin  (A61)  as 
shown  in  Figures  96  and  97. 


18 


iiiiiiiyiiiM 


SECTION  V 


ACOUSTIC  HELD  PREDICTION  METHOD 


5.1  htroduction 

A procedure  for  estimating  fiuctuating  pressures,  hereafter  referred  to  as  noise,  on 
STOL  aircraft  is  presented  in  this  section.  The  procedure  is  mainiy  concerned  with 
predictions  aft  of  the  nozzle  exit  plane  and  in  direct  view  of  the  engine  exhaust  flow 
stream,  in  region  A of  Figure  98.  An  approach  to  extending  the  procedure  to  indirect 
points,  in  Region  B,  is  aiso  provided. 

The  procedure  yields  1/3  octave  band  spectrum  estimates  associated  with  five  propul- 
sion/flap noise  sources  and  with  turbulent  boundary  layer  (TBL)  activity. 

The  total  noise  is  then  taken  to  be  the  (power)  sum  of  the  separate  source  spectra.  A 
typical  (low  speed/high  power)  situation  is  suggested  in  Figure  99. 

The  general  range  of  application  of  this  procedure  is  summarized  in  Figure  100,  and  is 
discussed  in  Sections  5.2-5.4. 

Methods  for  modifying  or  supplementing  the  1/3  octa/e  band  estimate  procedure  to  yield 
power  spectral  density  estimates  are  discussed  briefly  in  Section  5.5.  Finally,  comparisons 
between  measured  and  estimated  1/3  oct^e  band  noise  spectra  are  presented  in 
Section  5.6. 

Within  the  context  of  the  present  contract, (i)  the  exterior  surface  noise  is  taken  to  be 
the  principal  function  governing  airframe  vibration,  and  (ii)  exterior  surface  noise  is 
considered  to  be  known  no  better  than  is  the  airframe  vibration.  Exterior  surface  noise 
estimation  hence  becomes  a part  of  the  general  overall  problem  of  estimating  airframe 
structure  vibration  given  airplane  configuration  and  structure  details  and  given  engine 
and  airplane  operating  parameter  values.  Hence,  within  the  context  of  Phase  I of  the 
present  contract,  noise  estimation  is  considered  as  independent  of  the  vibration 
estimation  problem  in  which  the  noise  environment  Is  known,  and  vice  versa.  Develop- 
ment of  the  vibration  prediction  procedure  is  thus  broken  into  two  independent  parts,  (i) 
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vibration  response  prediction  given  the  noise  excitation,  and  (ii)  noise  excitation 
prediction  given  the  airpiane/engine/flap  configuration  and  operating  status.  The  focus  in 
this  section  is  on  the  iatter  prediction  procedure. 

The  procedure  presented  has  resulted  from  support  provided  both  under  the  present 
contract  eind  under  recently  completed  NASA  contract  NAS2-9328  (Reference  5). 
YC-14  data--which  was  the  prime  data  used  in  the  development  of  the  noise  estimation 
procedure  presented  herein~was  analyzed  and  a general  characterization  for  the  YC-  l^i 
airple^ne  developed.  This  characterization  was  then  generalized  and  formalized  into  a 
prediction  procedure  applicable  to  any  USB  STOL  airplane.  Without  question  the 
motivation  and  support  of  both  contracts  has  been  essential,  and  without  support 
provided  under  both,  the  development  of  the  procedure  presented  herein  would  not  have 
been  possible. 

3.2  Scope 

The  present  procedure  provides  1/3  octave  band  estimates  of  fluctuating  pressures  on 
USB  STOL  aircraft  surfaces  primarily  aft  of  the  nozzle  exit  plane,  and  in  direct  view  of 
(most  of)  the  engine  exhaust  field. 

The  estimate  for  a typical  field  point  P in  Figure  101  is  taken  to  depend  primarily  on  the 
characteristics  of  the  jet  flow  field  closest  to  P.  For  the  purposes  of  the  estimate,  a 
ribbon  Idealization  of  the  flow  field  is  employed.  The  procedure  is  specifically  oriented 
to  points  strongly  scrubbed  by  the  exhaust  flow  stream,  as  well  as  points  up  to  about  5 
(hydraulic)  nozzle  diameters  away  from  the  flow  boundary,  and  up  to  about  10  diameters 
downstream  of  the  nozzle.  The  procedure  is  most  applicable  for  cold  secondary,  dual 
flow  nozzles  (with  bypass  ratios  between  2 and  6 and  aspect  ratios  less  than  about  3) 
where  the  bottom  lip  is  integral  with  the  wing  top  surface. 

The  procedure  yields  1/3  octave  band  spectrum  estimates  for  each  of  the  following: 

o 3et  mixing  noise  in  the  presence  of  a scrubbed  wing/flap  system  with  or  without 
vortex  generators 
o Near-nozzle  noise 

o Trailing-edge  noise 

Noise  associated  with  (partial)  separation  of  the  exhaust  flow  from  flaps 


o 


o Turbulent  boundary  layer  noise 

o Exhaust  shock  noise 

A summary  flow  diagram  for  the  overall  procedure  Is  present  In  Figure  102.  Note  that 
turbomachinery  noise  (for  both  Inlet  and  exhaust)  Is  not  Included  In  this  estimation 
procedure. 

5 J Estimation  Procedure  Development 

As  noted  In  Section  5,2,  the  prime  source  of  data  drawn  upon  has  been  that  for  the 
YC*  and  which  is  summarized  and  discussed  in  Reference  5, 

The  most  obvious  characteristics  of  all  the  static  and  low-speed  YC-14  data  was  (a)  the 
simple,  single-peaked,  gently  rolling-off  spectrum  shape  of  the  noise  for  all  points  close 
to  or  scrubbed  by  the  jet  mixing  region  of  the  exhaust  flow  field,  and  (b)  the  inverse  ratio 
between  the  spectral  levels  and  their  distance  from  where  the  flow  field  roughly  seemed 
to  be.  The  general  shape  of  noise  spectra  is  illustrated  in  Figure  103.  The  dependence 
on  distance  away  from  the  flow  field  is  illustrated  In  the  same  figure,  and  also  In  Figure 
104  dc  103,  in  which  the  position  of  the  flow  field  (as  reflected  In  the  position  of  the  USB 
flaps)  is  changing.  The  effect  of  forward  velocity  Is  illustrated  In  Figure  106,  and 
suggests  a reduction  In  peak  spectral  level  and  an  Increase  In  the  frequency  of  the  peak 
level  with  Increasing  forward  velocity. 

On  the  basis  of  observations  as  these,  a flow  field  idealization  model  and  a jet  mixing 
noise  model  were  developed.  These  models  together  yielded  fuselage  field  point  to  flow 
boundary  separation  distances,  and  noise  levels  consistent  with  the  smoothed  behavior  of 
much  of  the  YC-14  ground  (and  some  low-speed  flight)  surface  noise  data. 

Estimate  procedures  for  the  remaining  noise  source  components  addressed  in  the  overall 
procedure  were  then  built  up  to  account  for  the  most  obvious  deviations  of  the  data  from 
the  estimated  jet  mixing  noise  component.  The  exhaust  shock  noise  component  was 
based  on  correction  of  distinctive  deviations  (re.  the  other  components)  observed  during 
high- speed/ high- altitude  (cruise)  operations. 


The  effort  summarized  in  the  above  three  paragraphs  is  discussed  in  detaii  in  Refer- 
ence 5.  In  particuiar,  extensive  discussions  of  the  characterization  of  the  fiow  field,  jet 
mixing  noise,  trailing  edge  noise,  exhaust  shock  noise  and  turbulent  boundary  layer  noise 
are  to  be  found  there. 

5A  Estimation  lYocedure 

Computational  aspects  of  the  estimation  procedure  are  divided  into  9 sections,  and  which 
are  presented  in  detail  in  Appendix  A of  Volume  II.  These  sections  arei 


Appendix  A Section 

Subject  Addressed 

3.1 

Characterization  of  the  Flow  Ribbon 

3.2 

Geometry  Computations 

3.3 

Jet  Mixing  Noise 

3.4 

Near-Nozzie  Noise 

3.3 

Trailing- Edge  Noise 

3.6 

Separation  Noise 

3.7 

Turbulent  Boundary  Layer  Noise 

3.8 

Exhaust  Shock  Noise 

3.9 

Estimation  for  Indirect  Field  Points 

Briefly,  the  parameters  used  to  characterize  the  flow  field  (ribbon)  idealization*,  shown 
in  Figure  107,  are  Its  maximum  width,  W*,  (or  angle,  8*,  and  its 

trail-off  angle,  8'.  These  are  computed  in  Section  3.1.  The  geometric  factors  which 
enter  into  these  computations  are: 

♦The  specific  idealization  is  as  follows:  The  flow  exits  the  nozzle  with  a width  equal  to 
the  nozzle  exit  width,  flush  with  the  wing  surface.  The  ribbon  spreads  linearly  in  width 
with  position  downstream  of  the  nozzle  exit  plane  until  it  reaches  the  beginning  of  the 
strongly  curved  portion  of  the  flap.  Thereafter  its  width  remains  constant,  and  its 
direction  of  flow  (as  viewed  from  above)  parallel  to  the  engine  centerline.  It  initially 
remains  attached  to  the  strongly  curved  portion  of  the  flap,  turning  to  the  angle  O',  at 
which  point  it  separates  from  the  flap  and  continues  on  a straight  course  at  the  elevation 
angle  O'. 
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0 


nozzle  side  lip  angles 
0 nozzle  top  lip  (kickdown)  and  bottom  angles 
0 wing  surface  inclination  angle 
0 nozzle  width,  height  and  effective  exit  area 
0 skew  angle  of  nozzle  exit  plane 

0 distance  from  nozzle  exit  plane  to  start  of  strongly  curved  portion  of  flap  system 

0 size  of  nozzle  side  door  opening  (If  present) 

In  addition,  the  following  operational  parameters  enter  Into  the  computationst 

0 static  flow  turning  angle  of  the  propulsion/flap  system  at  the  specific  flap  setting 
considered 
0 airplane  speed 

0 engine  exhaust  mixed  jet  velocity 

The  coordinates  of  the  field  point,  P,  at  which  the  noise  estimates  are  sought,  are  next 
computed  (Section  3.2)  in  terms  of  (see  Figure  101) 

6 s minimum  distance  of  P from  ribbon 
S « downstream  coordinate  of  P as  measured  along  ribbon 

Information  required  for  these  include  (In  addition  to  those  flow  ribbon  parameters  from 
Section  3.1)i 

0 Coordinates  of  the  field  point  P 

0 Coordinates  of  fixed  reference  point  P^  on  the  nozzle  exit  plane  (see  Figure  3 of 

Vol.  II,  Appendix  A) 

The  values  of  these  two  coordinates,  when  normalized  by  Dp|,  where 
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i 


i 

I 


and 


Agpp  a effective  area  of  engine  nozzle  exit  plane  (including  effect  of 
both  primary  and  fan  flows), 


along  with  the  values  of  parameters  listed  In  Figure  108,  are  then  used  to  compute  the 
estimates  of  the  various  noise  components) 

o jet  mixing  noise  (Section  3.3) 
o near-nozzle  noise  (Section  3.4) 

o trailing  edge  noise  (Section  3.3) 

0 separation  noise  (Section  3.6) 

0 exhaust  shock  noise  (Section  3.8) 

In  the  case  of  turbulent  boundary  layer  (TBL)  noise,  6 and  S are  not  required,  but  rather 
(see  Section  3.9) 

X a boundary  layer  growth  length  along  the  airframe  surface  to  the  field 
point. 

Operational  parameters  required  are 

V = representative  flow  velocity  along  boundary  layer  growth  path 

P a representative  flow  density  of  fluid  along  boundary  layer  growth  path, 

rather  than  those  given  in  Figure  108. 


Each  noise  component  is  characterized  in  terms  of 


0 a generalized  spectrum  shape 

0 spectrum  shape  peak  level, 

0 frequency  at  which  the  peak  level  occurs,  fp|^,  as  suggested  In  Figure  109.  In  the 
case  of  jet  mixing  noise  a modification  of  the  generalized  spectrum  shape  Is 
introduced  if  vortex  generators  are  deployed  into  the  flow  (see  Section  3.3). 
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All  generalized  spectrum  shapes  t . essentially  based  on  measured  1/3  octave  band  data. 
As  an  example  the  data  shown  previously  In  Figure  103  for  measurement  points  3,  7, 

13  and  14  was  in  part  used  to  define  the  generalized  spectrum  shape  for  jet  mixing  noise. 


Peak  level  computations  generally  combine  conceptual  scaling  rules  for  velocity  and 
density  effects  with  empirically  observed  effects  of  the  dimensionless  distance  of  the 
field  point  away  from  the  flow  ribbon,  6 /Dj^,  and  downstream  of  the  nozzle  exit  plane, 
S/D|^.  The  general  form  used  for  SPLpj^  Is 


In  which 


+ (6/Dpj)+A2  ^ 


P rep 


(p  I for  TBL  noise 

P|  I for  all  other  noise  components, 


and 


P 

Pj 

Po 


representative  fluid  density  for  TBL 
engine  mixed  exhaust  jet  density 
reference  density  at  which  SPL°p|^  is  defined. 


With  regard  to  the  second  term 


j for  jet  mixing  noise 


; for  near  nozzle  noise  and  separation  noise 

.4 


for  trailing  edge  noise 


for  TBL  noise 
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= 0 I for  shock  noise 


Vj  B engine  mixed  exhaust  jet  velocity 

s airplane  velocity 

= reference  velocity  at  which  SPLpj^  is  evaluated 

Cj  = engine  mixed  exhaust  sound  speed 

V = representative  velocity  for  TBL  activity 


The  third  and  fourth  terms  in  the  general  expression  for  SPL  > , i.e.,  Aj  and  are 
empirically  determined  relations  for  the  effect  of  and  S/Dj^  on  peak  level.  Recall 

that  is  the  engine  exhaust  nozzle  hydraulic  diameter.  Note  that  in  the  case  of  TBL, 
terms  of  the  form  of  Aj,  and  do  not  appear^  as  suggested  previously, 

Finally  the  fifth  term  in  the  general  expression  for  SPLpj^,  is  the  reference  peak 
spectrum  level  at  the  reference  conditions  of  s p^,  s V^,  and  at  6/D^»o,  and 
usually  at  S/Dj,^  = 3. 


The  computation  of  f^,.,  i.e.,  the  frequeiicy  at  which  the  peak  spectral  level  occurs,  Is 

Pi'  jh 

typically  a blond  of  empirical  relations  for  the  effect  of  '2nd  S/Dj^,  and  scaling 

rules  based  on  the  size  of  and  speed  at  which  eddies  closest  to  the  field  point  are 
generated  or  are  convected  past  the  field  point.  A general  form  is 


where 

r 


18V  /D 

I'  H } for  jet  mixing  and  shock  noise 
S/D^  + 3 


3.6 


for  near  nozzle  noise 
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U8  Vj/D|^  . jfQ,.  trailing  edge  noise 
+ 3 

J_  V|/6jj7}  for  separation  noise 


and  in  which  Vj,  S,  Dj^  and  V are  as  defined  previously,  and 


SjE  = 

6te  “ 
6tb  = 


distance  from  nozzle  exit  to  flap  trailing  edge  as  measured  along  the 
flow  ribbon 

distance  between  flow  ribbon  and  flap  trailing  edge 
thickness  of  TBL  at  field  point 


For  the  second  term  in  the  general  expression  for  fpj^ 


« 


for  jet  mixing, 
trailing  edge  arid 
shock  noise 


I 1 (for  all  other 

^ noise  components 

The  third  term  in  the  general  expression  for  fp|^,  namely  is  empirical  in  nature 

accounting  for  frequency  changes  (typically  Increasing)  associated  with  increasing 
distance  between  field  point  and  the  flow  ribbon  point  of  closest  approach. 


5.9  Application  to  Nan-ow  Band  Noise  Estimation 

The  procedure  referred  to  in  Section  3.4  and  in  Appendix  A of  Volume  II  was  developed 
using  1/3  octave  band  acoustic  data,  and  hence  was  itself  posed  in  1/3  octave  band 
terms.  However,  for  purposes  of  structural  vibration  estimation,  a power  spectral 
density  estimate  is  required  rather  than  a 1/3  octave  band  format. 
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Limited  examination  of  power  spectral  density  data  (corresponding  to  the  1/3  octave 
band  data  used  in  the  prediction  method  development)  shows  these  to  be  in  general 
smooth  curves  free  of  distinctive  (narrow  band)  peaks.  Hence  the  1/3  octave  band 
procedure  can  in  principal  be  quite  simply  extended  to  predict  exterior  surface  noise 
field  power  spectral  density: 

(a)  For  each  source  component  generalized  1/3  octave  band  spectrum  shape  curve 
denoted  as  spl(f),  determine  a generalized  power  spectral  density  shape  curve, 
denoted  by  psd(f),  such  that 

[ 

Kfj)  - I psd(f)df  = 0 

where 

i(f^)  s 

and  f|  is  the  i th  1/3  octave  band  center  frequency. 

(b)  determine  the  location  of  the  frequency  of  the  peak  of  the  generalized  psd  shape 
curve  with  respect  to  the  peak  of  the  generalized  spl  shape  curve. 

The  determination  of  a psd  function  yielding  a desired  set  of  l(f|)  values  at  f = f^ 
(i=l,2,...,n),  requires  one  to  assume  a form  for  the  psd  function  (containing  undetermined 
coefficients).  Values  for  ttriese  coefficients  are  then  chosen  to  make  the  right  hand  side 
of  the  equation  relating  l(f|)  and  psd(f)  as  close  to  zero  at  f3f|i  (lsl,2,...,n)  as  desired. 
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As  example,  a simple  form  for  a psd  function  is 


p8d(f)  = 


2 2 


2-1/6,  <,<2l/6f 


In  this  case  the  Cj's  are  given  simply  by 

1 

‘^i  = ^jT7r"pT7r^  fj 

This  format  Is  in  fact  used  in  Section  VI,  except  that  It  Is  applied  external  to  the 
estimate  procedure,  rather  than  within  it. 

^.6  Comparisons  of  Measured  vs  Estimated  Acoustic  Data 

Comparisons  are  presented  between  measurements  and  estimates  generated  with  the 
procedure  summarized  in  the  previous  section  and  presented  fully  in  Appendix  A of 
Volume  II.  Four  sets  of  YC- 1^  measurement  points/f light  conditions  are  considered: 

Set  1:  At  five  flap  measurement  points  all  at  the  same  STOL  approach  condition. 

Set  2:  At  nine  fuselage  measurement  points  at  the  same  brake  release  condition. 

Set  3:  During  various  phases  of  a take-off,  covering  brake  release  to  climbout,  at 

one  fuselage  location. 

Set  4:  At  various  extensions  of  the  USB  flaps  at  two  fuselage  locations. 

Values  of  geometry  parameters  for  the  YC-iii  used  In  arriving  at  these  estimates  are 
listed  in  Figure  110.  Measurement  point  locations  coordinates  are  summarized  in  Figure 
ill.  A general  diagram  showing  the  location  of  these  and  other  YC-14  measurement 
points  appears  In  Figure  40,  presented  previously  in  Section  3.6.  Values  of  operating 
parameters  for  the  10  ground  and  flight  conditions  examined  are  shown  in  Figure  112. 


29 


Finally,  a measurement- point/flight-condition  cross  reference  list  is  presented  in 
Figure  113. 

With  regard  to  each  of  the  resultant  estimates  presented  in  each  of  Figures  lU  - 139, 
the  measured  value  curve  Is  Indicated  by  solid  circular  symbols,  the  curve  of  (total) 
estimated  noise  by  solid  up-side-down  triangular  symbols.  The  remaining  operi  symbols 
indicate  values  of  various  noise  components,  as  noted  on  the  lower  part  of  each  figure. 
Note  also  that  the  lower  part  of  each  figure  provides  a brief  indication  of  the  operating 
status  of  the  airplane,  and  under  “NOTES**  the  location  of  the  measurement  point  and  the 
name  of  the  flight  condition.  With  regard  to  these,  the  following  abbreviations  are  used, 
ALT  = airplane  altitude 
SPEED  = airplane  speed 
N1  = engine  fan  shaft  rotational  speed 
VMIX  = engine  mixed  primary  and  fan  exhaust  velocity 
USBFA  = USB  flap  angle 

Under  Notes: 

BS  s body  station  location  of  measurement  point 
WL  = water  line  location  of  measurement  point 
BL  = butt  line  location  of  measurement  point 
VG  = vortex  generator 

The  above  selection  of  measurement  point/flight  conditions  spans  a reasonably  broad 
scope  of  the  STOL  airplane  low  speed  operations.  Note  that  this  selection  covers  an 
overall  acoustic  level  variation  of  about  30  dB  (i.e.,  128  to  159  dB)  and  a 10  to  i 
frequency  range  within  which  the  peak  spectral  level  falls  (i.e.,  UO  to  ^00  Hz).  Figures 
11^-127  give  an  indication  of  the  prediction  procedure  ability  to  assess  correctly  the 
effect  of  measurement  point  location  relative  to  the  engine  exhaust  stream  on  acoustic 
levels.  Figures  127-131  indicate  directly  the  procedure's  ability  to  handle 

forward  velocity  effects,  while  Figures  132-139,  USB  flap  position  effects. 
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Since  a portion  of  the  data  exhibited  in  these  figures  was  used  in  the  development  of 
! geneidlized  spectrum  shapes,  reference  levels,  etc.,  appearing  in  the  prediction  pro- 

‘ cedure,  the  comparisons  shown  Indicate  primarily  the  self  consistency  of  the  procedure. 

Based  on  comparisons  presented,  and  the  rtinge  of  locations  and  operating  conditions 
‘j  , covered,  it  is  felt  that  the  procedure  is  to  first  order  highly  self  consistent.  Improved 

self  consistency  could  be  achieved  with  further  evaluation  of  YC-l<i  and  other  (as  for 
example  QSRA)  data.  However,  such  an  effort  is  felt  to  be  beyond  the  scope  of  the 
present  program. 

Si" 
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PHASE  n PARAMETRIC  STUDIES 


Analytical  techniques  and  computer  programs  used  in  Phase  I were  extended  for  use  in 
Phase  II  of  the  program.  These  extensions  were  used  to  explore  the  variation  of 
environmental  vibration  2ind  acoustic  ieveis  in  STOL  type  aircraft  smaiier  and  larger 
than  the  medium  STOL  airplane  studied  in  Phase  1.  The  technicai  analysis  of  Phase  II 
consisted  of  three  parts:  (1)  deveiopment  and  eigenvalue  analysis  of  finite  element 
structural  models,  (2)  definition  of  the  acoustic  environment,  and  (3)  random  harmonic 
analyses.  The  USB  flap  structure  and  fuselage  structure  near  the  wing  root  were  the  two 
areas  studied  for  parametric  effects.  These  studies  foliow  in  Sections  6.0,  7.0  and  8.0. 


SECTION  VI 


FLAP  STRUCTURE  VIBRATION  PREDICTION- PARAMETRIC  STUDIES 
6.1  Small  STDL  - QSRA  USB  Flap 

The  QSRA  airplane  was  chosen  to  represent  the  small  scale  STOL.  Since  we  had  some 
experience  with  modification  of  the  airplane  for  NASA,  we  had  access  to  structural 
details  and  could  use  this  data  in  our  structural  representations. 

The  QSRA  USB  flap  structural  representation  was  derived  from  the  structure  shown  in 
Figures  1<»0,  141,  and  142.  The  upper  skin  material  was  .071"  AL  301. 

The  finite  element  models  of  the  QSRA  flap  were  developed,  as  shown  1n  Figure  143, 
with  31  node  points.  The  coordinates  for  the  nodes  were  assigned  as  shown  in  Figure  144. 
The  beam  elements  were  determined  as  shown  in  Figure  143  and  the  plate  elements  are 
given  in  Figure  146.  This  model  was  then  input  to  the  Structural  Analysis  Program  (SAP 
IV)  to  obtain  the  mode  shapes  and  frequencies  of  the  QSRA  flap.  It  is  to  be  noted  that 
the  response  of  the  six  locations  are  on  the  flap  structure  and  would  not  include  the 
attachment  point  (which  is  assumed  rigid  in  this  analysis  and  would  have  no  motion). 

The  first  20  frequencies  were  calcula+ed  and  are  listed  in  Figure  147.  The  mode  shapes 
are  shown  in  Figures  148  through  157. 

The  acoustic  input  was  determined  from  the  data  as  described  in  the  acoustic  parametric 
prediction  Section  Vlll  of  this  report.  The  excitation  points  of  the  flap  were  determined 
by  dividing  the  USB  flap  upper  surface  into  20  panels  as  shown  previously  in  Figure  146 
with  a given  pressure  acting  over  each  panel.  The  power  spectra  and  cross-power 
spectra  for  the  pressures  acting  on  the  panels  were  extrapolated  from  the  data  that  were 
given  for  the  six  locations.  The  acoustic  data  have  been  given  in  Section  VIII  in  the 
discussion  of  the  fluctuating  pressure  estimates  for  the  50,000  lb  STOL  airplane.  The 
extrapolation  of  this  data  for  all  the  panels  was  accomplished  using  the  extrapolation 
technique  used  for  the  YC-14  predictions.  The  frequency  responses  for  six  arbitrarily 
chosen  locations  on  the  flap  used  the  dynamic  analysis  computer  program  as  was  used  in 
the  YC-14  calculation  in  Phase  I.  The  locations  chosen  are  shown  in  Figure  158  as  points 
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1 through  6.  Structural  response  results  have  been  obtained  for  g s .06>  .09,  .12  and  .1^. 
Results  oi  g = .09  are  shown  in  the  plots  of  Figures  159  through  161  for  locations  1,  3, 
and  5,  for  the  airplane  condition,  STOL  approach,  Nj  = 85%,  50**  USB  flaps. 

A comparison  with  flight  data  for  the  QSRA  flap  is  given  in  Figure  162.  Location  1 was 
taken  as  the  point  for  comparison  to  the  accelerometer  that  was  mounted  on  the  flap 
actuator  (A13V). 

Lar  ge  STOL  > USB  Flap 

To  formulate  a design  for  the  large  STOL  airplane  flap  we  follow  the  scaling  chart  listed 
belows 


SMALL 

LARGE 

WT 

50,000  lbs 

1,000,000  lbs 

SPAN 

100  ft 

270  ft 

THRUST 

20,000  lbs 

400,000  lbs 

The  linear  scale  factor  of  large  STOL  airplane  to  small  STOL  airplane  is  2.7.  If  we 
scaled  the  QSRA  flap  accordingly,  the  large  STOL  airplane  flap  would  have  a chord  of 
130  inches  and  span  of  190  inches. 

The  design  of  the  large  STOL  flap  would  differ  only  slightly  from  the  QSRA  for  purposes 
of  this  study.  The  details  of  design  are  shown  in  Figures  163,  and  16U  with  the  values  of 
the  components  listed  in  Figure  164.  Figure  165  indicates  the  6 locations  chosen  for  the 
analyses  solutions. 

The  large  STOL  USB  flap  model  was  input  to  the  SAP  IV  program  with  node  points  as 
shown  in  Figure  166,  the  coordinates  as  shown  in  Figure  167  the  plate  elements  as  shown 
in  Figure  168  and  the  structural  component  values  as  shown  in  Figure  169. 

The  output  from  the  SAP  IV  program  is  listed  in  Figure  170  for  the  first  20  Modes.  The 
mode  sl-iapes  are  shown  in  the  following  Figures  171  through  177. 
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The  acoustic  input  was  obtained  as  detailed  in  Section  8.4.  The  response  plots  of  the 
large  STOL  USB  flap  model  to  the  STOL  condition  of  50®  flap  setting  and  85%  power  are 
shown  in  Figures  178  through  189. 

We  sec  the  response  of  the  large  USB  flap  has  lower  frequency  content  than  the  QSRA 
flap  with  the  response  of  the  64  Hz  mode  and  the  second  and  third  modes  clearly  seen  in 


SECTION  VII 


i FUSELAGE  STRUCTURE  VIBRATION  PREDICTION-PARAMETRIC  STUDIES 

I 

I 7.1  SmaU  STOL  - QSRA  Fuselage 

V ! 

I The  QSRA  fuselage  model  was  chosen  to  represent  the  upper  fuselage  from  the  top 

surface  of  the  wing  to  airplane  center  line.  The  area  chosen  can  be  seen  in  Figure  190 
! and  191.  Several  models  were  constructed,  one  being  a 16  node  model,  where  response 

calculations  were  made.  However,  for  better  accuracy,  a larger  model  was  finally 
I selected  as  shown  In  Figure  192  that  had  77  nodes.  The  coordinates  for  this  model  were 

selected  as  given  in  Figure  193  with  the  beam  and  plate  elements  given  in  Figure  194  and 
t 195.  The  actual  stringer  used  in  the  QSRA  airplane  are  shown  in  drawings  of  Figure  196 

I i and  197  which  were  Included  ir  'he  model  calculations.  Figure  198  summarizes  the 

i'  structural  values  used  in  the  QSRA  fuselage  model. 

i: 

F' 

The  results  of  the  calculation  of  SAP  IV  with  the  Berm- Plate  QSRA  fuselage  model  of  70 
f;i  nodes  Is  listed  in  Figure  199  for  the  first  20  modes.  The  shapes  of  the  lower  modes  have 

f'v  been  drawn  in  Figures  200  through  209. 


The  model  was  then  excited  by  the  acoustic  Input  described  in  Section  8.3  for  the 
locations  shown  and  for  a damping  value  of  .09.  The  results  are  shown  in  Figures  210, 
211  and  212. 

The  fuselage  model  represented  the  upper  fuselage  structure  from  airplane  structure, 
body  station  345  to  450  and  WL  198  to  209.  The  QSRA  accelerometers  that  would 
correspond  to  this  region  are; 

A5V,  A6L  BS  400  Side  Frame/Strlnger  Junction 

A7V,  A8L  BS  500  Ceiling  Longeron 

A comparison  with  flight  data  for  the  QSRA  fuselage  is  given  in  Figure  213  where 
location  7 of  the  math  model  Is  compared  to  accelerometer  A66  which  was  mounted  in 
the  QSRA  airplane  at  BS  400  on  a side  frame-stringer  junction.  The  levels  show 
satisfactory  agreement.  Location  1 was  also  compared  to  test  data  and  indicates  the 
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test  data  to  be  lower  than  predicted  up  to  approximately  70  Hz.  (Figure  214).  This  area 
Is  somewhat  out  of  the  direct  Impingement  area  and  could  be  responsible  lor  some  error 
in  the  assumed  excitation  or  the  model. 

7JL  Large  STOL  - PuMiage 

The  fuselage  model  was  scaled  up  from  the  QSRA  data  with  length,  area  and  mass 
factors  proportional  to  2.3,  (2.5)^,  and  (2.5)^  The  scaled  structural  values  were  input  to 
the  SAP  IV  program  for  modal  frequencies  listed  In  Figure  213  and  plotted  for  mode 
shapes  shown  In  Figure  216  through  221. 

The  acoustic  Input  was  then  determined  as  given  In  Section  VDI  and  the  structural 
response  for  the  stringer  frame  locations  are  given  In  Figures  222,  223  and  224. 
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NOISE  FIELD  PARAMETRIC  PREOICTiON 

This  section  discusses  fuselage  and  flap  surface  noise  predictions  which  were  generated 
to  provide  an  excitation  input  source  to  the  vibration  analyses  of  "smalF*  and  "large" 
STOL  airplanes  discussed  in  Sections  VI  and  VII.  The  selected  airplane  geometries, 
operating  conditions,  field  point  locations  and  necessary  scale  factor  relations  for 
dimensions  and  operating  parameter  values  are  set  forth.  An  example  noise  prediction 
tabulation  and  plot  for  the  small  airplane  are  presented  and  discussed.  Their  relation  to 
the  prediction  procedure  of  Section  Vl  is  discussed.  The  simple  manner  In  which  the 
estimates  for  the  small  airplane  can  be  applied  to  the  large  airplane  are  stated. 

8.1  STOL  Airplane  Prediction  Parameters 

The  basic  airplane  geometries  for  which  flap  and  fuselage  surface  noise  levels  have  been 
generated  (and  which  have  previously  been  discussed  in  Sections  VI  and  VII  aret 

o A four-engine,  50,000  ib  gross  weight  airplane  - the  QSRA  configuration/design  is 
used. 

0 A four-engine,  1,000,000  ib  gross  weight  airplane  - a scaled  version  of  the  QSRA  is 
used  for  simplicity. 


The  two  operating  conditions  for  each  are  chosen  asj 


Brake  release  (100%  rated  thrust) 

Airplane  speed  (V^)  = 0 

Engine  mixed  jet  velocity  (V^)  = 870  ft/sec 

USB  flaps  at  0°  (fully  retracted) 

STOL  operation  (85%  rated  thrust) 

Airplane  speed  = 110  ft/sec  (65  knots) 
Airplane  altitude  (ALT)  = 6500  ft 
Engine  mixed  jet  velocity  = 680  ft/sec 


USB  flaps  extended  to  50 
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Fixed  vortex  generators  are  assumed  for  both  airplanes  at  both  operating  conditions, 
following  the  scheme  actually  used  for  QSRA. 

Figures  22?  — 230  in  the  present  section  describe  the  50,000  lb  QSRA  type  airplane, 
and  the  field  points  at  which  estimates  have  been  made.  Geometric  data  shown  in  these 
figures  is  based  on  QSRA  drawings,  primariiy  Boeing  Dwg.  340-000003. 

The  i, 000 ,000  ib  gross  weight  airplane  was  taken  to  be  a scaled  up  version  of  the  above 
50,000  lb  airplane,  with  its  dimensions  proportional  to  the  cube  root  of  the  ratio  of  the 
gross  weights.  Hence  all  dimensions  of  the  QSRA  airplane  apply  to  the  1,000,000  lb 
airplane  upon  multiplying  by  a scale  factor  SF,  of 

5F  n ^^7^0,000/50,000  s 2.71 

A second  assumption  made  is  that  the  engines  of  the  1,000,000  lb  airplane  are  exact 
scaled  replicas  of  those  of  the  50,000  lb  airplane,  and  all  havp  the  same  engine  cycle, 
by-pass  ratio,  etc. 

Under  these  assumptions  it  Is  further  assumed  that  engine  mixed  jet  velocity,  airplane 
speed  and  USB  flap  angle  are  the  same  for  both  airplanes  at  brake  release  at  100%  rated 
thrust,  and  at  STOL  operation  at  85%  rated  thrust. 

8.2  Rpediction  Method 

The  method  described  in  Section  V was  used  to  generate  the  estimates  for  the  small 
STOL  airplane,  as  well  as  providing  the  simple  guidelines  needed  to  apply  these  to  the 
large  STOL  airplane. 

Through  Boeing  in-house  support,  a computerized  version  of  the  current  procedure  has 
been  developed  for  the  CDC  systdm.  This  program  USBESTO),  generates  tabulated 
spectral  value  lists  by  noise  component  and  in  totai,  as  weli  as  piotting  files.  Via  an 
existing  computer  plotting  program  constructed  as  a part  of  a 1976-1977  AFFDL/NASA 
contract  effort  to  measure  YC-i4  cabin  noise,  graphs  of  the  estimates  can  also  be 
generated.  These  programs  have  been  used  to  generate  all  estimates  appearing  in  this 
report. 
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An  example  of  a computer  tabulation  and  corresponding  computer  plot  are  shown  in 
Figures  231  and  232«  respectively. 

The  one  principle  difference  between  the  computerized  version  of  the  noise  prediction 
procedure  and  that  presented  in  Appendix  A of  Volume  II  is  the  definition  of  zones. 
Within  the  program,  and  on  the  output  tabulation  forms  (e.g.,  Figure  231)  the  following 
definitions  are  applied: 

Zone  0:  All  points  above  the  upper  wing  surface  which  are  forward  of  the  engine 
nozzle  exit  plane. 

Zone  1:  All  points  above  or  on  the  wing  upper  surface  which  have  S values  between  O 

and  1#^. 

Zone  2:  All  points  above  or  on  the  wing  upper  surface  which  have  S values  between 
and  S*. 

Zone  3:  All  points  above  or  on  the  wing  upper  surface  with  an  S value  greater  than  S', 

and  all  points  with  an  S value  greater  than 

Zone  4:  All  points  v/ith  an  S value  less  than  and  which  are  below  the  wing  lower 

surface. 

The  definitions  of  S,  S'  and  S^.^  used  above  are  the  same  as  those  used  in  Appendix  A of 
Volume  II. 

Symbols  and  abbreviations  appearing  on  the  tabulation  and  plot  forms  of  Figures  231  and 
232  generally  follow  or  are  mnemonics  fcr  those  used  in  Appendix  A of  Volume  II.  For 
instance,  with  regard  to  Figure  231: 

ALT  = airplane  altitude 
VA  = airplane  forward  speed 

V3  = engine  mixed  (primary  and  fan)  exhaust  flow  velocity 
VGS  s vortex  generators 

R/RO  s ratio  of  at  altitude  air  density  to  sea  level  air  density 
THETAS  = 0j,  per  Appendix  A of  Volume  II 
THETAP  =9',  per  Appendix  A of  Volume  II 
EX  = "exit" 

TE  = "trailing  edge" 

DELTA  = 6,  per  Appendix  A of  Volume  II 
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Abbreviations  for  noise  components  include: 

MIX  = jet  mixing  noise 

NN  = near  nozzle  noise 

TE  = trailing  edge  noise 

SEP  = separation  noise 

TBL  = turbulent  boundary  layer  noise 

SUM  = noise  associated  with  power  sum  of  all  above  noise  components 

5.3  Fluctuating  Ptesaure  Estimates  - Small  STOL 

A complete  set  of  tabulations  and  plots  for  the  8 flap  and  wing  field  points  and  the  9 
body  field  points  are  included  as  Appendix  B of  Volume  II.  A summary  of  overall  levels 
at  each  field  point  at  each  of  the  two  conditions  considered  (brake  release  and  STOL 
operation)  is  shown  in  Figure  233. 

Note  in  this  figure  that  separate  estimates  due  to  the  inboard  pair  of  engines  alone  and 
due  to  the  outboard  pair  of  engines  alone  are  included;  The  estimate  procedure  is 
designed  for,  and  is  based  on,  data  from  airplanes  with  two  symmetrically  placed 
engines,  as  the  YC-i4.  Hence  treatment  of  a four-engine  airplane  must  be  handled 
indirectly,  i.e.,  two-engines  at  a time.  To  this  end  it  is  assumed  that  for  noise  purposes 
the  contributions  from  each  pair  of  engines  can  be  treated  independently,  and  then 
summed  on  a power  basis  to  obtain  total  noise.  (Measured  QSRA  does  exist  for 
evaluating  this  assumption  at  least  for  noise  on  the  fuselage,  but  such  a check  has  not 
been  made,  See  for  instance  Boeing  Document  D6-^f7I18,  "QSP.A  Flight  Test-Noise," 
J.  E.  Sommers  and  A,  J.  Bohn,  14  Dec.  1978.)  The  values  in  Figure  233  and  in  the 
Appendix  tabulations  and  plots  (see  Volume  II)  indicate  that  under  this  assumption,  the 
contribution  of  the  outboard  engine  pair  never  contributes  more  to  the  total  (summed) 
noise  on  the  fuselage  or  inboard  flaps  than  3 dB  and  in  most  cases  less  than  1 dB, 

8.4  Fluctuating  Pressure  Estimates  - Large  STOL 

Because  of  the  very  special  relationships  imposed  between  the  two  "paper"  STOL 
airplanes  which  are  considered,  i.e.. 


o They  are  geometrically  similar  with  each  dimension  of  the  large  airplane  being  2.71 
times  the  corresponding  dimension  of  the  smaii  eurnlane. 

o The  operating  parameter  vaiues  are  exactly  the  same  for  both  at  brake  release,  and 
then  again  at  the  STOL  operation  considered, 

the  spectral  estimates  for  the  small  airplane  can  be  extended  with  but  small  error  to  the 
large  airplane.  Specifically, 

o overall  levels'  and  spectral  levels  for  large  airplane  s overall  and  spectral  levels  for 
small  airplane 

0 Frequencies  (and  frequency  scales)  for  large  airplane  = (dimensional  scale  factor)"  ^ 
X frequencies  (and  frequency  scales)  for  small  airplane. 

That  is,  the  spectra  for  the  large  airplane  are  those  for  the  small  airplane  upon  dividing 
the  frequency  scale  (of  the  small  airplane  spectra)  by  the  dimensional  scale  factor,  being 
2.71  for  the  two  airplanes  considered  here. 

These  simple  spectral  relationships  arise  from  the  interaction  of  the  similar  geometries 
and  identical  operating  condition  values  for  the  two  airplanes  considered  with  the 
following  properties  of  the  noise  estimation  method  summarized  In  Section  V and 
described  in  detail  in  Appendix  A of  Volume  11. 

0 All  component  noise  levels  depend  upon  the  ratio  of  field  point  distance  to  nozzle 
hydraulic  diameter,  which  for  the  airplanes  under  study  are  the  same  for  both  (see 
in  particular  Section  5.4). 

0 Beyond  this  ail  levels  depend  in  addition  on  airplane  and  engine  operating 
parameter  values  (i.e.,  airplane  altitude  and  speed,  engine  mixed  exhaust  velocity, 
li  sound  speed,  and  density,  and  USB  flap  angle),  and  which  for  the  airplanes  under 

study  are  the  same  for  both  (see  in  particular  Section  5.4). 
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o With  the  exception  of  TBL  noise,  characteristic  frequencies  scale  inversely  with 
engine  hydraulic  diameter,  which  for  these  two  airplanes  are  related  by  the  scale 
factor  of  2.71. 

In  the  case  of  TBL  noise,  the  characteristic  frequencies  scale  inversely  to  boundary  layer 

O 

thickness  which  for  the  two  airplanes  considered  here  go  as  (geometric  scale  factor)  = 

e 

(2.71)  = 2.2.  At  brake  release  wher^  the  airplane  speed  is  essentially  0,  this  effect  is 

negligible.  At  the  STOL  operatibn  condition  the  TBL  noise  component  is  small  compared 
to  the  jet  mixing  noise  component  so  that  again  the  non-simple  scaling  effect  for  TBL 
has  a negligible  effect  on  the  total  noise.  (See  in  particular  Section  A.4.7  of  Appendix  A 
of  Volume  II.) 

The  manner  chosen  for  relating  the  geometries  and  operating  statuses  of  different 
weight  airplanes  in  this  study  leads  to  probably  the  most  concise  relationship  possible 
between  surface  noise  fields  for  different  airplanes.  To  first  order  the  scaling 
approaches  employed  seem  quite  reasonable.  Hence,  to  first  order,  surface  noise  levels 
would  be  expected  to  remain  about  the  same  for  airplanes  of  different  size,  (but  with  the. 
same  number  of  engines  of  the  same  by-pass  ratio),  while  the  characteristic  frequency  of 
the  noise  would  become  lower  with  increasing  airplane  size. 

A detailed  examination  of  differences  in  missions,  aerodynamic  and  propulsive  perform- 
ance of  airplanes  of  differing  weight  would  undoubtably  lead  to  less  simple  geometric 
and  operations  relations.  In  turn,  these  would  lead  to  less  easily  descrlbable  differences 
in  the  surface  noise  fields.  However,  within  the  limitations  of  the  noise  prediction 
procedure  which  has  been  developed  in  the  present  study,  such  differences  in  geometry 
and  operations  should  be  addressable  directly,  and  without  difficulty. 
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SECTION  IX 


ENVIRONMENTAL  VIBRATION  PREDICTION  COMPARISONS 

Comparison  of  the  applicable  military  standards  was  made  here  to  Indicate  the  impact  of 
the  predicted  vibration  levels  In  STOL  aircraft.  Pour  pages  of  MIL-STD  SIO  C showing 
the  predicted  levels  for  this  type  aircraft  are  shown  In  Figures  234  through  237.  Figure 
238  upper  curve  Indicates  the  levels  in  MIL-STD  810  C for  1 HR  test  and  the  lower 
curve,  the  actual  environment  based  upon  a 20,000  hour  life  as  specified  in  test  factor 
from  Reference  4,  Shock  and  Vibration  Handbook,  p.  2k-2k, 

The  plot  of  the  STOL  responses,  Figure  239,  compared  to  the  MIL-STD  810  C levels  show 
the  lower  frequencies  of  the  large  STOL  airplane  to  be  of  some  concern  since  this  energy 
can  be  transmitted  into  primary  structure  at  these  frequencies.  The  levels  shown  for  the 
small  STOL  in  Figure  240  will  have  considerable  attenuation  as  we  move  from  the 
strlnget  frame  locations  down  to  heavier  frame  support  structure  where  equipment  would 
be  located,  but  MIL-STD  810  C still  would  appear  to  be  somewhat  inadequate  for  the 
smaller  STOL  aircraft  In  the  frequency  range  from  125  to  300  Hertz. 
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Figure  2.  YC-14  USB  Flap  Simulation  for  the  EKS-SAP  IV 


Figure  4.  USB  Flap  SAP  IV  Data  Card  laage 
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MODE 
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(RAO/SEC) 
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(CYCLES/SEC) 

PERIOD 
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1 

1.9S11E*02 

3.1689E401 

3.1557E-02 

2 

2.3100E'»02 

3.87B5E401 

2.7200E-02 

3 

2.5561C402 

A .0682C401 

2.R561E-02 

A 

R.9966E-»02 

7.9S23E*401 

1 .2575C-02 

5 

7. 7028E+02 

1.2259E402 

&.1570E-03 

6 

8*047Q£-»02 

1.28C7E402 

7.fi081E-03 

1 

1. 06A<tE«03 

1.&9A0E402 

5.9032E>03 

B 

1.A597E+03 

2.3232E402 

^.30^'»E"03 

9 

1.6A65E403 

2,6205E402 

3.6160E-03 

10 

1.7853E403 

2.8098E402 

3.6592E-03 

Figure  5.  USB  Flap  Mode  Frequency  Spectrum 
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Figure  7.  USB  Flap  Mode  Shapes 


I I I I I I I I • * 

U^t^tJWJUWU^Li 
sBair)«dii'i(v<rrf«'a 
WiriVlMAlidNlNff^*' 
v^^<i4e4BM<^<or*M 
rjn««Mr%r*^Na«« 
n«4<4iA^c<^in#4r 

I I II 


i«»i*oooaoo 

I I 4 I I I I I 
iUlWU  k.iUWUiUJ 

‘<MN.#n«a*»H(rN 

III  1 


l#V1PlK1W*4 


I I I I I 1 I I I I 
U w WW  UU  LI  u W 

f%virgiN9«c<MK«n 

r)inin*^r4ii30rf)f*««i 


>00000  OOOOI 


I I I I I I I I I I 
UlUJ^UlU  UwUMU 


I I I I.  I I I I I I 
UUU«JWULIUUUI 


OOoooooeoa 
I I I I I I .1  I I I 
UUWLiuUlUUJiJU 
• «rooi>iinr*»-»<r 

III  II  II 


*5  N >0  ftf  «W  N to  N <M  IS 

to  oooaoooooo 


MMiCirinNP«nur«4>  ^MV«Uifs9^0«'^>«i4 

d9fSp49«#sorsvi*«ri 

p>*r«M^>or«Ma^ci  ^noa»r>«i^Oin 

‘!sc^in<soo»«tnMh»  f^oo'csOMO<o*<4e^ 


I I I I I I t I I I 

y kiuiLiyuyiyyy 

^KtrONirlOOO^ 

•^yr«fs<A«p<ionr* 


I I I I II 


oowr«oa*«o*H  o 

I I * I I • I t I * 

yLiywuwsjuiwu 

OrflO^rflK^wl-^vO 
a 0inu*to*«vfitaW)o 

.or>>n*^yis<oo*>««4 
• I I 


ra  iirt  «'  ts  ^ I 

II  II 


I ^ wi*  a a 

> o o cf  a o o a 

I I I « I I ♦ 

I u y w y ui  ui  w 

‘ QH  m ^ 

I o N 4 o « rs 

I a fv  fs  1*1  i-t  ^ 

’ a Cl  •«  4 cr  o ^ 

I Cl  a «*>  **  <0  v« 


(StAOfSOv  fSI^Oc^ 
II  I II  1 


-a  a^4J 

V Dies 
X C*r- 
•r-*r-  O 
U.XO. 


IS  10  M «•  V 4M  a 
II  I 11 


I I 1 I « I I I 

uuutyuuiyu 

cswiNrOi^Sl^a 

moiot^atoo^fs 

rr^ooOp^i/ia 

t II  II 


Figure  7.  USB  Flap  Node  Shapes  (Continued) 
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Figare  ?.  USB  Flap  Mode  Shapes  (Continued) 


Figure  7.  USB  Flap  Node  Shapes  (Continued) 
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Figure  7.  USB  Flap  Mode  Shapes  (Cont.) 
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USB  Flap  Hode  Shapes  (Cont.) 


Figure  7.  USB  Flap  Mode  Shapes  (Concludad) 


Figure  9.  USB  Flap  Mode  2.  36.765  Hz 
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Figure  11.  USB  Rap  Node  6,  128.07  Hz 


Figure  12.  IBB  Rap  Hade  8.  232.32  Hz 
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Figurt  15.  Map  of  Upper  Triangle  of  the  CPSO  Matrix  for  Model  I 
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Figure  16.  Exaaple  of  Print-Out  for  Three  Accel erooeter  Locations  for  Model  I 


DIAGONAL  CPSD  MATRIX 
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(10-3) 

a^/RAD  aVHz 

2 
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(10-3) 

G^/RAD  G^/HZ 

3 

141 

(10. 

GVRAD 

28 

GVH2 

26 

.044 

.28 

.145 

.91 

.088 

.055 

27 

.059 

.37 

.194 

1.22 

.130 

.82 

28 

.061 

.51 

.267 

1.68 

.205 

1.29 

29 

.114 

.72 

.360 

2.39 

.363 

2.28 

30 

.170 

1.07 

.588 

3.69 

.821 

5.16 

31 

.322 

2.02 

1.247 

7.84 

3.410 

21.43 

32 

.636 

4.00 

2.604 

16.36 

8.171 

51.34 

33 

.660 

4.15 

2.234 

14.04 

2.172 

13.65 

34 

1.177 

7.40 

3.836 

24.10 

1.966 

12.35 

35 

2.821 

17.72 

9.104 

57.20 

3.725 

23.40 

36 

11.265 

70.78 

36.255 

227.80 

13.920 

87.46 

37 

29.75 

186.92 

96.70 

601.30 

37.15 

233.42 

38 

6.53 

41.03 

21.06 

132.32 

8.98 

56.42 

39 

3.27 

20.55 

10.61 

66.66 

5.74 

36.07 

40 

5.51 

34.62 

18.11 

113.79 

12.68 

79.67 

42 

2.74 

17.22 

9.02 

56.67 

5.82 

36.57 

43 

1.34 

a-.Az 

4.41 

27.71 

2.56 

16.08 

45 

.648 

4.07 

2.13 

13.38 

1.07 

6.72 

50 

.307 

1.93 

1.009 

6.34 

.417 

2.62 

60 

.192 

1.21 

.620 

3.90 

.236 

1.48 

70 

.433 

2.72 

1.260 

7.92 

.323 

2.03 

75 

1.960 

12.32 

5.38 

33.80 

.898 

5.64 

76 

3.217 

20.21 

8.783 

55.19 

1.356 

8.52 

77 

5.97 

37.51 

16.21 

101.85 

2.345 

14.23 

78 

13.048 

81.98 

35.230 

221.36 

4.867 

30.58 

79 

30.497 

191.62 

82.315 

517.20 

11.002 

69.13 

80 

33.669 

211.55 

90.780 

570. 

11.978 

75.26 

81 

16.47 

103.48 

44.44 

279. 

5.88 

36.95 

lAeter  Response  Predictions*  G".03 

i9 


FREQUENCY 

(HZ) 

DIAGONAL  CPSD  MATRIX  | 

1 

1421 

(10“*) 

6^/RAO  G*/HZ 

— 

2 

1417 

6VRA0  G^/HZ 

3 

1428 

(io“3) 

g^/rao  g^/hz 

85 

2.455 

15.43 

6.686 

42. 

.994 

6.25 

90 

.939 

5.90 

2.610 

16.4 

.346 

2.17 

100 

.454 

2.85 

1.310 

8.23 

.281 

1.77 

no 

.329 

2.07 

.984 

6.18 

19.35 

121.58 

120 

.350 

2.20 

1.676 

10.53 

.143 

0.90 

122 

.500 

3.14 

3.42 

21.49 

.141 

0.89 

123 

.591 

3.71 

3.98 

25.01 

.144 

0.90 

124 

.604 

3.80 

3.48 

21.87 

.149 

0.94 

125 

.622 

3.91 

3.101 

19.48 

.160 

1.01 

126 

.697 

4.38 

3.223 

20.25 

.182 

1.14 

127 

.823 

5.17 

3.774 

23.7 

.220 

1.38 

128 

.859 

5.40 

4.020 

25.26 

.246 

1.55 

129 

.660 

4.15 

3.080 

19.35 

.210 

1.32 

130 

.459 

2.88 

2.034 

12.78 

.159 

1.00 

135 

.270 

1.70 

.833 

5.23 

.085 

.53 

140 

.279 

1.75 

.718 

4.51 

.076 

.48 

150 

.424 

2.66 

.805 

5.06 

.068 

.43 

160 

1.387 

8.71 

1.974 

12.40 

.068 

.43 

165 

4.967 

31.21 

6.554 

41.18 

.087 

.55 

167 

10.53 

66.16 

13.73 

86.27 

.122 

,77 

168 

15.52 

97.52 

20.206 

126.96 

.154 

.97 

170 

20.05 

125.98 

26.18 

165. 

.188 

1.18 

171 

15.50 

97.39 

20.32 

128. 

.161 

1.01 

175 

3.995 

25.10 

5.401 

34. 

.087 

.55 

185 

0.758 

4.76 

1.148 

7.2 

.066 

.41 

195 

0.358 

2.25 

.598 

3.8 

.065 

.41 

Figurfi  17.  YC-14  Flap  Accelerometer  Response  Predictions,  G-.03  (Concluded) 
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Figure  26.  Natural  Nodes  of  Hedlua  STOL  Flaps,  Model  II 


Figure  27.  YC-14  USB  Flap  Finite  Elenent  Model  II  Mode  1 Frequency  = 31.88  Hz 
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Figure  35.  Node  9 Frequency  > 159.63  Hz 


Figure  36.  YC-14  USB  Hap  Finite  Element  Model  II.  Mode  10.  Frequency  = 210.49  Hz 


Figure  37.  Mode  11,  Frequency  = 275.09  Hz 


Figure  38.  Node  12  . Frequency  - 296.80  Hz 


Figure  39.  Ibde  13  , Frequency  - 313.25  Hz 
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Figure  43.  USB  Flap  Acoustic  Excitation  Spectrui,  Microphone  N4Q 
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Figure  46.  YC-14  USB  Flap  Response  Predictions 
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Figure  58.  Low  Frequency  Fuselage  Model  1 
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Figure  65.  Fuselage  Skin  Model  II 
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Figure  68.  Natural  Modes  of  Low  Frequency  YC-14  Fuselage  Model  I 
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Figurt  77.  Nicrophont  D«t«  for  Nodlun  STOL 
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Calculated  Response  of  Hodels  I and  II  for  Stringer  Acceleroaeter  ASS 
Coapared  to  Flight  Test  Data 


Accel erowter  No.  1481 


Calculated  Response  of  Models  1 and  II  for  Fuselage  Skin  Accelcroaetcr  A61 
Coapared  to  Flight  Test  Data 
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NUMBER 
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( RAO/SEC } 
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PERIOD 
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1 
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Figure  81.  USB  Flap  Model  II  Frequency  Spectrum 
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Figure  83.  Fuselage  Model  II  Modal  Plots,  Frequency  > 145.57  Hz 


Figure  85.  Fuselage  Model  II  Modal  Plots,  Frequency  * 234.64  Hz 


Figure  86.  Fuselage  Model  II  Nodal  Plots,  Frequency  ■ 284.06  Hz 
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Figure  88.  Fuselage  Model  II  Modal  Plots,  Frequency  « 341.17  Hz 
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Figure  91.  High  Frequency  Model  III  Frequency  Spectrum 
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Figure  93.  Fuselage  Model  III  Modal  Plot,  Frequency  « 326.01  Hz 


Figure  94.  Fuselage  Model  III  Modal  Plot,  Frequency  - 376.41  Hz 


ee  D3  Accelerometer  No.  1474 


Figure  96.  USB  Flap  Model  III  Response  Comparison  With  Stringer 
Accelerometer  {A58)  for  Three  Sets  of  Assumed  Damping 
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USB  Flap  Model  III  Response  Conparlson  with  Skin  Mounted 
Accelerometer  (A61}  for  Three  Sets  of  Assuaed  Oaaping 


Figure  98.  General  Regions  for  Application  of  USB  STOL  Airplane  Estimation  Procedure 
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Figure  99.  Gener«1  Arrangement  of  Component  Noise  Source  Estimates  Making 
Up  Total  Noise  Estimates  (For  Typical  Cllmbout  Condition) 


• APPLIES  TO  ANY  USB/STOU  AIRPLANE  CONFIGURATION  WITH 

• NOZZLE  FLUSH  TO  WING  SURFACE 

• LOW-SPEED  OPERATION  FROM  TAKEOFF  THROUGH  CRUISE  (VaA/J  < 1) 

• NOZZLE  ASPECT  RATIO  UP  TO  ~ 8 

• BYPASS  RATIO  FROM  ~ 2 TO  -6 

• ACCOUNTS  FOR 

• ENGINE  MIXED  JET  VELOCITY  (Vj) 

• ENGINE  MIXED  JET  DENSITY  (p{) 

• AIRPLANE  FORWARD  SPEED  (Va) 

• SIZE  AND  POSITION  OF  NOZZLE  WITH  RESPECT  TO  FUSELAGE 

• NOZZLE  CONFIGURATION 

• US8  FLAP  ANGLE 

• WING/FLAP  CONFIGURATION 

• VORTEX  GENERATORS  (IF  PRESENT) 

• NOZZLE  SIDE  DOOR  (IF  OPEN) 

• FLOW  TURNING  CAPABILtTY  OF  NOZZLE/FLAP  SYSTEM 


Figure  100.  Scope  of  USB/STOL  Aircraft  Noise  Estimation  Procedure 
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Figure  101.  Conceptual  Relation  of  Typical  Field  Point  P to  Flow  Field 
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Figure  102.  Summary  Flow  Diagram  for  USB/STOL  Aircraft 
Fluctuating  Pressure  Estimation  Procedure 
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Figure  104.  Effect  of  USB  Flap  Position  on  Exterior  Fuselage 
Overall  Levels 
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Figure  10b.  Effect  of  USB  Flap  Position  on  Exterior  Fuselage  Spectra 


- Hirplane  forward  velocity 

Vj  engine  mixed  exhaust  jet  velocity 

6py  » static  flow  turning  capability  of  propulsion/flap  system 
(when  trailing  edge  flap  system  is  at  6f;degrees) 

pj  » engine  mixed  exhaust  jet  density 

Cj  a engine  mixed  exhaust  sound  speed 


Figure  108.  Airplane  Operating  Parameters  Used  In  Estimation  Procedure 
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Figure  110.  YC-14  Geometry  Parameter  Values  Used  In  Noise  Prediction 
Procedure  Exercise 
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Measurement  Locations  Used  for  Prediction  Procedure  Exercise 


63 


CONDITION  NUMBER 


Figure  113.  Neasurement  Point/Operattng  Condition  Cross  Reference  List  for 
Prediction  Procedure  Exercise 
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Figure  118.  Results  for  M41  at  Condition  3160 
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Figure  121.  Results  for  MOB  at  Condition  7132 
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Figure  125.  Results  for  M15  at  Condition  7132 
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Figure  130.  Results  for  H20  at  Condition  7135 


2 I49III0'  2 94SII  10^  2 9 4 S • I 10 

riCOUCNCY  !HZI 


PLOT 

iimoL 

f 

y 

□ 

0 

A 

X 


Z-OUCCR  CONO.  AIT.  SPEED  NI  VHI X USBFA  OVERALL  ' 

NO.  - HQ.  I F T I I FPS  l LfiEtU  LfPlJ  ULfii  .IPBI  ■ 


M20 

7I3B  100  220  3710  lOSO  0 

143 

M20 

7131 

143 

M20 

7131 

IIS 

M20 

7I3B 

0 

M20 

7131 

IIS 

M20 

7131 

130 

M20 

7130 

143 

EXT  loor  ism  vlzsa  blozi 

PREOICICO  TOTAL  NOISE  .CREATED 
PREOICICO  TBL  NdiSI 
PREOICTEO  SEP  NOISE 
PREOICICO  COCE  NOISE 
PREOICICO  NN  NOISE 
PRIOICIEO  MIXING  NOISE 


CLIMB. 13$  KNTS 
7S/0//I0. 
79/07/10. 
71/07/10. 
79/07/10. 
79/07/10. 
79/Or/IO. 


Ftgur«  131.  Rtsults  for  N20  tt  Condition  7)36 


12 


MICRO- 


PREO  Cl  ON  PROCEDURE  DEMONSTRAl  ON 


1 '4  S • • tO^  2 3 4 

6 16 

lo’ 

2 3 

4 s e • 10 

rRCOUCNCY 

IHZI 

X-OUCER  CONO.  AIT.  SPEED 

Nl 

VMIX 

USBFA 

OVERALL 

- NO.  NO.  .IfTI.  LLPJJ 

IRPMI 

].rp.s,.i 

ijim 

-IDii- 

MI3  7196  10000  :;2S 

29S0 

8 30 

9 

147 

MI3  7196 

148 

MI3  7196 

129 

M13  7196 

0 

MI3  7196 

102 

M13  7196 

134 

M13  7196 

148 

CXr  PAIR  •S67S  VL22D  BL102 

FLAP 

CYCLE. VC 

OVN 

PREOICTEO  tOIAL  NOISE .CREATED 

79/07/10. 

PREOICTEO  I9L  NOISE 

79/07/10. 

PREOICTEO  SEP  NOISE 

79/07/10. 

PREDICTED  EDGE  NOISE 

79/07/10. 

PREDICTED  NN  NOISE 

79/07/10. 

PREDICTED  HIXINC  NOISE 

79/07/10. 

Figure  132.  Results  for  M13 

it  Condition 

7196 

PREDICTION  PROCEDURE  DEMONSTRATION 


V W W l-J 

34SCII0^  2 34 

5(9 

10^  2 3 4 S • • 10 

PRCOUCNCY 

IHZI 

X-DUCCR  COND.  ALT.  SPEED 

HI 

VMIX  USBFA  OVERALL 

HO.  NO.  LFJI  irPS.1 

LRPHJ 

ILP.S.I  UUJU  ...1 DLI  - 

H13  7193  11000  210 

2950 

830  29  144 

H13  7193 

145 

K13  7193 

112 

M13  7193 

102 

H13  7193 

115 

MI3  7193 

no 

HI3  7193 

145 

EXT  FAIR  BS87S  VL220  9LI02 

FLAP  CYCLE. VO  OWN 

PREOICTEO  TOTAL  NOISE  .CREATCO 

79707/10. 

PREOICTEO  I6L  NOISE 

79/07/10. 

PREOICTEO  SEP  NOISE 

79/07/10. 

PREOICTEO  EDGE  NOISE 

79/07/10. 

PREOICTEO  NN  NOISE 

79/07/10. 

PREOICTEO  MIXING  NOISE 

79/07/10. 

Figure  133.  Results  for  N13  at  Condition  7193 

PREDICTION  PROCEDURE  DEMONSIRA T ION 


2 

W kJ  IvJ  A M M W 

s^seiio’  2 34S$I 

10*  2 3 

4 5 6 8 10^ 

PitCOUCNCV  IT4ZI 

PLOT 

X-DUCCR  COND  ALT.  SPCCD  HI 

VMIX  USBFA 

OVERALL 

ium 

-AO... _->LQ...  IfTI.  1.IPJJ  iifMJ 

UAii  iJEC..I 

JJllL- 

• 

HI3  7I9S  10000  215  2950 

830  41 

144 

T 

Mt3  7I9S 

143 

□ 

M13  7195 

112 

0 

MI3  7195 

114 

0 

HI3  7195 

116 

A 

Ht3  7195 

126 

X 

MI3  7195 

143 

No^rs 

EXT  FAIR  BSB75.  V1220  BL102 

FLAP  CYCLE. VG 

UP 

T 

PRCbiCTCO  TOTAL  NOISE  .CRCATCO 

79707/10. 

□ 

PRCOICTEO  TBL  NOISE 

79/07/10. 

0 

PREDICTED  SEP  NOISE 

79/07/10. 

0 

PREDICTED  EDGE  NOISE 

79/07/10. 

A 

PREDICTED  NN  NOISE 

79/07/10. 

X 

PREDICTED  MIXING  NOISE 

79/07710. 

Figure  134.  ' Results  for  M13  at  Condition  7195 

185 

PREDICTION  PROCEDURE  DEMONSTRATION 


Sii: 


i 9 

4 s e 

• 10'  234 

SIS 

IC*  2 9 

4 S • 1 10 

FREOUENCY 

IH2I 

PLOT 

X-OUCCR 

COND.  ALT.  SPEED 

HI 

VMIX  USBFA 

OVERALL 

SYMBOL 

--  WO. 

.....NQ.  irn  LtfiiJ 

IRPMl 

JJ/-S.1  UULIU 

f 

M13 

7192  10700  213 

2950 

930  70 

140 

? 

HI  3 

7192 

143 

P 

Ml  3 

7192 

112 

5 

Ml  3 

7192 

135 

0 

MI3 

7192 

lOB 

A 

Ml  3 

7192 

124 

X 

M|] 

7192 

142 

WQ^CS 

CXT  FAIR  BSI7S  VL720  BII02 

FLAP  CYCLE. VG 

UP 

? 

PRCOICTCO  TOUL  NOISE  .CREATED 

79/07/10. 

□ 

PREOICTCO  TBl  NOISE 

79/07/10. 

$ 

PREOICTEO  SEP  NOISE 

79/07/10. 

0 

PREOICICD  EDGE  NOISE 

79/07/10. 

A 

PREOICIED  NN  NOISE 

71/07/10. 

X 

PREOICTED  mixing  NOISE 

79/or/io. 

Figurt  135. 

Results  for  M13  et  Condition  7192 

OCTAVC  BANS  f ' SB  BE 


PREOlCnON  PROCCOURE  DCMUNSIRADON 


2 

94SIII0'  2 34 

SIR 

lO-’  2 3 

4 5 I • 10- 

mCOUCNCY 

IH2I 

i 

PLOT 

X-OUCCR  CONO.  alt.  SPCCO 

Nl 

vMix  user A 

OVERALL 

uim 

NO.  HO.  TTTI  liiU 

IRPMI 

IFP51  1 DEC  1 

1 DR  1 

f 

Mi4  7193  HOOD  210 

2990 

130  29 

149 

▼ 

MIA  7193 

147 

(P 

M14  7193 

112 

MIA  7193 

91 

0 

MIA  7193 

100 

A 

MIA  7193 

131 

X 

MIA  7193 

141 

lUl|Li 

(XT  lOOr  BS97S  VLIBO  BLI07 

FLAP  CYCLE. VG 

OWN 

? 

PRCOICTCO  total  NOISE .CftCATtO 

79/07/10. 

Q 

PRCOICTCO  TIL  NOISE 

79/07/10. 

f 

PNCOICTEO  SEP  NOISE 

79/07/10. 

0 

PRCOICTEO  EOGE  NOISE 

79/07/10. 

A 

PREOICTEO  NN  NOISE 

79/07/10. 

X 

PREDICTED  MIXINC  NOISE 

79/07/10. 

Figure  137.  Results  for  N14  «t  Condition  7193 

1/3  octave  IAND  Ot  »E  ^ 5 


PRE0IC110N  PROCCOURE  OEMONSIRAIION 


4 9 6 


2 3 4 9 1 

rncoucNCY  ihzi 


4 9 6 


PLOT  X-OUCER  CONO.  ALT.  SPEED  N1  VMIX  USBFA  OVERALL 


MttftOi, 

? 


NO.  NQ._  JJJJ,  Lf£Li  LREiil  LfP-SJ  IPTCI  1 D9 1 - 


nss  toooo  2IS 

7I9S 

7195 

7195 

7195 

7195 

7195 


930  41 


Kayi 

T 

□ 


EXT  OODV  BS675  VLISO  BLI07 
PREDICTED  TOTAL  NOISE  .CREATED 
PREOICTEO  TBL  NOISE 
PREOICTEO  SEP  NOISE 
PREOICTEO  EOCE  NOISE 
PRCOICTED  NN  NOISE 
PREOICTEO  MIXINC  NOISE 


FLAP  CYCLE. VC  UP 
79707/10. 
79/07/10. 
79/07/10. 
79/07/10. 
79/07/10. 
79/07/10. 


Figure  138,  Results  for  M14  at  Condition  7195 


/IT* 


DC  TATE  lARD  lH 


PREDICTION  PROCEDURE  DEHONSTRAT ION 


4 S I I IO< 


4 5 1 ■10'' 


ricouCNcy  ihzi 


HOT 

sthbol 

f 

T 

□ 


X-OUCCR 


CONO.  AIT.  SPEED  Nl  VMI X USBEA 

HC.  TETI  lEPSl  1,RP>H  ltJ>Sl  IJECJ 

7IB2  10700  213  2950  830  70 

7192 
7192 
7192 
7192 
7182 
7192 


458  8 10 


OVERALL 


T 

□ 

t 

A 

X 


EXT  BODY  IS87S  VllBO  OLIO. 
PREOICIED  TOTAL  NOISE .CREATED 
PREDICTED  T8l  NOISE 
PREDICTED  SEP  NOISE 
PREDICTED  EDGE  NOISE 
PREDICTED  NN  NOISE 
PREOICIED  MIXING  NOISE 


FLAP  CYCLE. VC  UP 
79/07/10. 
79/07/10. 
79/07/10. 
79/07/10, 
79/07/10. 
79/07/10. 


Figure  139.  Results  for  N14  «t  Condition  7192 


qSRA  USB  FLAP 


SCALE  1"  ■ 8 


22 

STR.  UPPER 

.071  X 

2,70  X 17,2 

AL  301  1/2  HARD 

26 

STR.  UPPER 

.071  X 

5.20  X 5.80 

AL  301  1/2  HARD 

27 

STR.  LOWER 

,071  X 

3,70  X 17.2 

AL  2024  -0 

28 

STGR  LOWER 

.071  X 

2.70  X 18,0 

AL  2024  -0 

31 

STQR  LOWER 

.071  X 

5,20  X 18,0 

AL  2024  -0 

4 

SKIN  LWR 

.071  X 

52,0  X 76,0 

AL  2024  ■ T3  SHEET 

3 

SKIN  UPPER  AFT  .071  x 

12,0  X 76,0 

AL  301  1/2  HARD 

2 

SKIN  UPPER  CTP 

t ,071  X 

32,5  X 76,0 

AL  301  1/2  HARD 

32 

SPAR  FRONT 

.071  X 

11,3  X 76,0 

AL  301  1/2  HARD 

33 

SPAR  REAR 

.071  X 

8.60  X 76.0 

AL  301  1/2  HARD 

F1gur«  141.  QSRA  USB  Flap  Schematic 


COORDINATES 


<2>  FORWARD 

AX 


16 

15 

8 

14 

.9  10 

13 

uUL 

12 

11 

10 

9 

< 

13  14 

15  ^ 

17 

18 

ORIGIN 

8 

7 

6 

5 

20  2W 

22 

23  24 



4 

3 

2 

1 

M 

28 

-29 



® ATTACHMENT  POINT  (PINNED) 


5 +25 

12  +16 


26  -16.0 
31  -25.0 


NODE  1 2 3 4 5 6 7 8 9 10  11  12  13  14 


X +25 


+25  k6.C 


iil 


Y l37.(i-18.q  0 H8.q37.0h37.0H21.5k8.0  0 18  21.5  37.0  +37. 


NODE  15  16  17  18  19  20  21  22  23  24  25  | 26  | 27  | 28  | 29 

X 0 0 0 0 0 -16  

Y 1-180]  0 -18. ( -21.5-37 03>37.0^21.5H8.0  0 ■la.(-21. 

Figure  144.  QSRA  USB  Flip  Finite  Elenents 

19S 


(LUM-22,  27) 
FSCAR 

(-32)  -26  -31 


3 1 3 


STRINGER 


(-26,  -31) 


DESCRIPTION 


LE 

4 

FS 

r 

STRINGER 

RS 

] 

TE 

1.5194 

0.3834 

1.1360 

0.0026 

3.0726 

10.9755 

0.8023 

0.4615 

0.3408 

0.0888 

10.0992 

0.6682 

1.4768 

0i5680 

0.9088 

1.2423 

0.2022 

0.0811 

0.6106 

0.2698 

0.3408 

.001027 

1.5717 

1.0428 

1.5194 

0.3834 

1 

1.1360 

0.0028 

1 

2.9358 

7.9337 

Figure  145.  QSRA  USB  Flap  Element  Properties 


1 


196 


11 


12 


V 

h' 


I 


I 


I 


17 


Figure  146. 


Figure  148.  Saill  USB  Flap  Model.  Modal  Plot  Frequency  ■ 150.44  Hz. 


Figure  152.  Snil  USB  Flap,  Hedal  Piet,  Frequency  ■ 306.38  Hz. 


Figure  155.  Small  USB  Flap,  Modal 
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165.  Large  STOL  USB  Rap 


Hsdcl  Node  Ntots 


USB 


PART 

-3 

SKIN, UPPER 

.250  X 140  X 190 

-4 

SKIN,  LOWER 

.250  X 130  X 190 

-22 

STRINGER,  UPPER  (Z) 

.1875  X 4.0  X 190 

■•■23 

STRINGER,  NOSE  (L) 

.1875  X 2.0  X 190 

-26 

STRINGER,  UPPER  (HAT) 

.1875  X 6.0  X 190 

-27 

STRINGER,  LOWER  (Z) 

.1875  X 4.0  X 190 

-31 

STRINGER,  LOWER  (HAT) 

.1875  X 6.0  X 190 

-32 

SPAR,  FRONT 

.375  X 14.0  X 190 

-33 

SPAR,  REAR 

.375  X 10.0  X 190 

Figurt  169.  Structural  Conponents  of  Largo  Airplane  USB  Flap 
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MODE 

CIRCULAR • 
FREQUENCY 

FREQUENCY 

PERIOD 

NUMBER 

(RAO/SEC) 

(CYCLES/SEC) 

(SEC) 

1 

4.02T«Ee02 

b.4097Ee0l 

•oisbo 

2 

5.b2»9E«02 

2.9491E431 

.01117 

3 

b.29aoe^o2 

1*0024E902 

.00999 

4 

7.24I1E«02 

l,lS2!ii€«02 

•009b9 

S 

«««b70E4-02 

l.3476Ee02 

.00742 

b 

U0b22E403 

4.bU05Ee02 

.00592 

? 

t.201SE4-0S 

1 .OIZUEOOB 

.00523 

h 

1..224Tt:’»'0i 

1.9q92Ef02 

.00513 

4 

1 ,2IS9E4*03 

2.D274E402 

,00493 

n 

URlSTE-rOJ 

2.2499E4^02 

.00444 

n 

1.5b27e't’03 

2.47}2E^C? 

•00405 

12 

uaeTTEt-os 

2,90BBfc>02 

,00.344 

13 

U^557C+03 

3.1  125E402 

.00321 

14 

2^0bSbE+03 

3.3l97E<-02 

,00301 

15 

2,3l.94E'«-Oi 

3.bVl5E402 

, 00271 

lb 

2,503SE403 

4.1U7E402 

.00243 

1? 

2,b5SAE<»03 

«.22b2E-f02 

.00237 

15 

2,7l25E'f03 

4.3l70Ef02 

.00232 

19 

2.9559e+03 

4,7013E^02 

.00213 

20 

3«2l29E+03 

5.1134E«02 

.00196 

Figure  170.  Print  of  Frequencies  for  Large  Airplane  USB  Flap  Model. 
3S  Nodes.  No  Camber 
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171.  Large,  §TDL  USB  Flap  Nodal  Plot,  Frequency  » 64.097  Hz. 


Figure  172.  Urge  STOL  USB  Rap  Itedal  Plot.  Frequency  . 89.43  Hz. 


Figure  173.  Large  STOL  USB  Flap  Nodal  Plot,  Frequency  » 100.24  Hz. 


Figure  175,  Largs  STOL  USB  Rap  Modal  Plot,  Frequency  * 134.76  Hz. 


Figure  176.  Large  STOL  USB  Flap  Modal  Plot,  Frequency  « 169.05  Hz. 


FREQUENCY  - H2(X10  ^ ) 

Figure  17B.  Response  Prediction  of  Large  STOL  USB  Flap  Position  2 
Damping  G-.06 


2H/C0S  5) 


2H/tDS 


335 


0 


50 


100 


2S0 


300 


Figure  185. 


150  200 

FREOLENCY  - HZCXIO  ‘ 1 

Response  Prediction  of  Large  STOL  USB  Flap  Position  6, 
Damping  G-.09 
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1 


RUTO  PSD  - £&  SQ)/RZ 


C&  SQ)/HZ 


1 

USB  LARGE  FLAP 

35  NODES 

G - .09 

POSITION  a 

Figure  193.  QSRA  Fuselage  Structure  Model  Coordinates 


qSRA 


.080  RAD. 


0.599 

0.581 


ACTUAL  SIZE 

I 


SCALE:  TWICE  SIZE 
MATERIAL:  AL  ALLOY  2024-T4 
SPECIFICATION:  qQ-A-267  TEMPER  T4 


Figure  196.  Skin  Stringers  for  QSRA  Fusel ege  Model 
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SCALE:  TWICE  SIZE 
MATERIAL:  AL.  ALLOY  2024-T4 
SPECIFICATION:  QQ-A-267  TEMPER  T4 


Figure  197.  Skin  Stringer  for  QSRA  Fuselage  Model 
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MODE 

NUMBER 

CIRCULAR 

FREQUENCY 

(RAD/SEC) 

FREQUENCY 

(CYCLES/SEC) 

PERIOD 

(SEC) 

TOLERANCE 

1 

l,0579F+03 

1 .fablBE^O? 

•0OS05 

6.9166e>15 

2 

1 ,7352E+03 

2.7bl6E*02 

.00302 

9,69HbF-15 

3 

1 .79C'ic*l*3 

2.Ha9l*E*02 

.110351 

n,65n£E-15 

4 

2.l00ir>03 

3.3“3^E+02 

,00299 

n. 

5 

2.16l!je+u3 

3,47?<iE+02 

.002BB 

1 .6763^-14 

h 

2,1925EVn3 

3. 4 49 4 E+ 02 

.(TD2B7 

1.8599F-14 

7 

2.2b5AEH-03 

3,5902E+02 

,00279 

1,7S70E-14 

6 

2.5214E+C3 

4.0137E+02 

,00249 

9.3718E-15 

9 

7.A2H1I4-03 

4.162TE4fi2 

,00239 

1 ,2945F*r4 

10 

2.7in9E+03 

4,5145E402 

.00232 

3,6496F-14 

1 1 

2.774‘>t+0i 

4,41S9E*02 

.00226 

1 ,1bl4F-14 

12 

^.ff9u2EtG3 

a.6nh?E+ft2 

,00217 

7,U5BE-15 

13 

3.4452t'^<-i 

S.4635F.+02 

.00162 

1 .bOSbE"*!  4 

14 

4,?25bC+0S 

«».7251E+02 

,00149 

5.1976E-10 

rs 

4, 4347 E *05 

7..0581E+n2 

.0-0142 

5. V59BE-09 

i»* 

4.51C0E+03 

7.1779E+02 

,00139 

1 ,60SflF-09 

17 

4,R202E>03 

7,353?E+02 

.00136 

1 .?OhlF-07 

l« 

4.6772F4^03 

7.4O39E^02 

.00134 

2.2859t-0h 

19 

4.t.999E  + &3 

7.4601E+02 

.00134 

l,0639K-10 

20 

4.7799r+03 

7,6074E+02 

.00131 

n,0063F-07 

\SPPU  Pnnnos  on  nctwvjiut  ccustfks 

,n679365150396E+07  . JO«nB7«239529SE^07  ,3236042419«590E’*'07 

,64235B7a04933lE't>07  ,fc97b63Sfl376031  Et07  .742249497  1 ?fl89E+07 

.lf»03319fa442l74E+0B  . 19B6336)  525223E+0B  .20543641  357BblE-^OR 


«r\/ 


Figure  204.  QSRA  Fuselage  Modal  Plot 
Mode  5 

f = 347.24  Hz. 


Figure  208.  QSRA  Fuselage  Modal  Plot 
Mode  9 

f » 418.27  Hz. 


flUTO  PSD  - t&  S0)/HZ 


ZH/tOS  0) 


flUTO  PSD  - (&  SO) /HZ 


ZH/I'DS 


QSRB  flight  2 

aSL  THKECFF 

RMS  3.4703 
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qSRA  FUSELAGE  PREDICTION 
Location  7 (Side  of  Model) 
G - .09  Assumed  Damping 


qSRA  BS400  Side  Frame 
Stringer  Junction 


Figure  213. 


FREQUENCY  - HZ(X10  J ) 

DELTA  F = 2.4111' 

Cotuparlson  of  QSRA  Fuselage  Response  with  Predicted  Values 
Location  7 


MODE 

lUNBER 

1 

CIRCULAR 

FREQUENCIES 

(RAD/SEC) 

3,b872ef02 

FREQUENCY 

(CYCLES/SEC) 

S.8bB3E401 

PERIOD 

(SEC) 

.01704 

TOLERANCE 

6,8504E-15 

2 

6.69131402 

1.0649E402 

.<'0939 

8.3P046-.15 

3 

7,4?92E+02 

1,1904E402 

.00340 

6.659bE-l 5 

4 

9,09u7E+02 

1.4466E4D2 

.00691 

9.0156E"15 

5 

9,17i»4E4U2 

1.4602E402 

.006SS 

2.2130E<«14 

b 

9.696'4£402 

1,5«35E402 

.00648 

0. 

7 

9,9966E402 

l,590»E402 

.00629 

7.4S72E"IS 

8 

1 ,OS7BC40S 

l,6»26E40a 

.00594 

6.b662E>*iS 

9 

1 ,I7r5E403 

l.bTdJ  E*02 

.00534 

1 ,0747fe-l4 

10 

1 .SOb'lE+o: 

2.397CE4C2 

.np'4l7 

4.9599E-12 

11 

1 ,629AE403 

2,5933E*02 

.003rt6 

1 ,1 IbOEtl 0 

12 

I .oagne+os 

2.61  ST £4  02 

.^0362 

2,!29aE-n 

13 

1,6H46E403 

2.ttttllE402 

.00373 

1 ,2010F-10 

i,7i'49y4e4u3 

2,7S42E4(»2 

.n03S9 

1,1S20E-.09 

15 

1 .«^i7-4t>tj3 

2,«76T>e^t>2 

.00345 

4.l99bF-I0 

lb 

1,95<i3E403 

3.110<iE4O2 

.no3»i 

l,693bE-09 

17 

2,047bL403 

3,2587E402 

.00307 

6,6344E'<07 

1» 

2,2l3aL403 

3,5221E40E 

.00264 

1 ,I237F«07 

19 

2.2b«6E403 

3,6042E402 

.00277 

3,59()2E-06 

20 

2.319QC403 

3,6906E4U2 

.00271 

2,U97E^06 

JPPtH  UttU'IOn  ElCfcNVALUt  CL'tSTf^S 

slSTSlOaseSfraiTE  + Ob  .«5a207«30  3B«i0aE+Oti 
,Uafl8<i<iaa')e»fe71E  + 07  .1«00«37  47  55';2PE>07 
,3a9')E7»»ai‘J9*77<JEt07  . 3B57fc7 1 .131 9r)38E  + 07 


.56‘i90a293<l<>7aOE  + Ofe 
,22ci0938l7fl0ql3E  + 07 
»a?3405felS#.?fr55Et07 


Figure  215.  Print  of  Frequencies  for  Large  STOL  Airplane  Fuselage  Model 
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Mode  2 f = 106.49  Hz 


Figure  Z19.  Large  STOL  Fuselage  Modal  Plot 
Mode  4 f = 144.68  Hz. 


Figure  220.  Large  STOL  Fuselage  Hodal  Plot 
Hode  5 f = 146.02  Hz. 


ISO  ailO  250  200 

QUEMCY  - HZCdO  M 

ITOL  Airplane  Fuselage,  Position  1 

’3 


Cl>  SQ)/HZ 


FREQUEMCr  - H2(X10  ' ) 

Figure  224.  PSD  Plot  of  Large  STOL  Airplane  Fuselage,  Position  7 
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rplane  - Top  View 


Figure  226.  QSRA  Type  Airplane  - Side  View 


View-UIng/Flap  GeMMtry  (Inboard  Engine) 
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HZ 
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116. 
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1£5. 

C»e; 

^<>.4  ■ 

131.3 
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1£9. 
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1£5. 
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£000. 

l£8. 
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1 £4 . 

84. 
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OO  k 

l£8.3 
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1£5. 

11£. 
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1££. 
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1£6.9 
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1£3, 
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1£1. 
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6*3. 

1£0. 
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1£3.9 

ORSPL 
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144.5 
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Figure  231 „ Computer  Prediction  of  Noise  for  QSRA  Type  STOL  Airplane. 
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Figure  232. 
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PRCOICTCO 


TOTAL  NOISE. CRCATCO 
TBL  NOISE 
SEP  NOISE 
EDGE  NOISC 
NN  NOISC 
MIXING  NOISC 


Prediction  for  QSRA  Type  Airplane. 
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79/03/21. 

79/03/21. 

79/03/21. 

79/03/21. 

79/03/21. 

79/03/21. 

USB"50  - Inboard  Engine 


, ACCELERATION  POWER  SPECTRAL 

REDUCTION  IN  DENSITY  IG*/HJt 


RANDOM  VIBRATION  ENVELOPE 


REDUCTION  FACTOR  FOR  MASS  LOADING 

Figure  237.  Mll-Std  810C  Random  Test  Levels 


287 


ACCELERATION  POWER  SPECTRAL  DENSITY  ^ 


MIL  STD  810 C 


/ 


m 


/\ 

/ \ 


g«.09  Assumed  Structural  Damping 


Large  STOL  Fuselage  (From  Fig  222) 


L—  Large  STOL  Flap  (From  Fig  179) 


[l)Th1s  curve  represents  the  MIL  STD  810C  level  extrapolated  to  the 
actual  environment  of  20.000  hr.  airplane  as  per  test  duration  fiictor 
of  shock  and  vibration  handbook.  Harris  & Credo  Page  24-24. 


FREOLENCY  - HZ 

Figure  239.  Comparison  of  Large  STOL  Response  to  Predicted  MIL-STD  610C 
Environmental  Vibration  Levels 
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FREOUEhICY  - HZ 

Figure  240.  Coinperlson  of  Small  STOL  Response  to  Predlctid  MIL-STD  810C 
Environmental  Vibration  Levels 
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